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Abstract 

In this paper we study surface operators in A/" = 4 supersymmetric Yang-Mills theory. 
We compute surface operator observables, such as the expectation value of surface opera- 
tors, the correlation function of surface operators with local operators, and the correlation 
function of surface operators with Wilson and 't Hooft loops. The calculations are per- 
formed using three different realizations of surface operators, corresponding respectively to 
the gauge theory path integral definition, the probe brane description in AdS^ x and the 
"bubbling" supergravity description of surface operators. We find remarkable agreement 
between the different calculations performed using the three different realizations. 
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1. Introduction and Summary 

In gauge theories, nonlocal operators can be understood as operators that insert a 
probe into the theory. Well known are the Wilson and 't Hooft operators, which are 
characterized geometrically by a curve in spacetime. A Wilson loop inserts into the theory 
a point-like charged particle while an 't Hooft operator inserts a magnetic monopole. These 
operators play an important role as order parameters of various phases of gauge theory. 

A surface operator on the other hand is characterized geometrically by a surface E in 
spacetime. This surface corresponds to the worldsheet spanned by a probe string that the 
surface operator inserts in spacetime. A surface operator 0-£ measures the response of the 
external string with which the field theory is probed, and therefore these operators may 
enlarge the list of known order parameters in gauge theories. 

In A/" = 4 SYM two types of super symmetric surface operators have been studied. One 
is akin to the 't Hooft operator in that it is an operator of disorder type, characterized by a 
certain supersymmetric vortex-like singularity near a surface E . One can also construct 
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a supersymmetric surface operator of order type , akin to a Wilson loop operator, which 
inserts into the A/" = 4 SYM path integral the WZW action integrated over E together 
with an associated Chern-Simons term, where the WZW action is constructed out of the 
A/" = 4 SYM gauge field along E (see also the work in 0)0. 

In this paper we go beyond identifying and classifying these operators and perform 
three classes of computations with the surface operators 0-£. We first study properties 
associated to the surface operator Ox; itself, namely evaluating its expectation value. Then 
we consider the interaction of a surface operator Oy: with local operators, including chiral 
primary operators and the stress-energy tensor. And then we analyze the interaction 
between a surface operator Oj: and other non-local operators in the gauge theory — 
certain supersymmetric Wilson loops and 't Hooft loops. These correlators exhibit a rich 
dependence on the parameters that a surface operator 0-£ is characterized by, which we 
review in Section 2.1. 

We study the disorder surface operators of [|l[] by using three alternative realiza- 
tions of them. In the first one we use the gauge theory definition in terms of supersymmetric 
codimension two singularities in the A/" = 4 SYM path integral In the second, a surface 
operator Oy: is realized in terms of a configuration of D3-branes ending on the boundary of 
AdS^ X along the surface E, which realizes the boundary conditions corresponding to a 
surface operator Oj: in the gauge theory living on the boundary. The third description of 
surface operators O-^ is in terms of "bubbling" geometries — ten dimensional solutions of 
Type IIB supergravity which are asymptotically AdS^ x S^. They were proposed in ^ as 
the gravitational description of surface operators and are related to the solutions of Lin, 
Lunin and Maldacena (LLM) ^ by analytic continuation. 

While these correlation functions are computed using a variety of techniques, ranging 
from purely field theoretical to string theoretic calculations using strings, branes and su- 
pergravity solutions, the results are very coherent. A rather complete picture arises that 
allows us to study the properties of surface operators in A/" = 4 SYM in different regimes 
of the gauge theory parameter space. It is quite remarkable that the results of our compu- 
tations computed in various regimes of A/" = 4 SYM show agreement with each other. The 
computations performed using semiclassical techniques in the gauge theory agree with the 
computations we perform using a bulk dual string theory computation. 

^ Disorder surface operators with higher pole singularities have been constructed in Q while 
rigid disorder surface operators have been studied in Q. Disorder surface operators have also 
been considered in Q. 
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The three reahzations of a surface operator O-^ are vahd in different regimes of the 
gauge theory parameter space. The semiclassical gauge theory description is appropriate at 
weak 't Hooft couphng A <^ 1. The probe D-brane description of a surface operator is vahd 
when ^ 1, A ^ 1 and when the number of D3-branes making up the corresponding 
surface operator is smaU. FinaUy, the supergravity description in terms of "bubbhng" 
geometries is appropriate when ^ 1, A ^ 1 and when the surface operator 0-£ is 
characterized in the probe description by a large number of D3-branes, in which case 
the backreaction of the D3 branes cannot be neglected. Nevertheless, surprisingly, we 
obtain quantitative agreement when the the same computation involving surface operators 
is performed using the different realizations of 

Summary of Results 

Before getting into the details of the calculations in the following sections, let us 
elaborate a bit on the quantities we compute in this paper. 

The first quantity we calculate is the expectation value of the surface operator 0-£. 
In this paper we consider maximally supersymmetric surface operators O^,, which can 
have the geometry of the plane or of a sphere, that is T, = or S'^. Naively, one may 
think that in a conformal field theory like A/" = 4 SYM, that the expectation value of 
Oj: cannot depend on the radius a of the 5"^ on which 0-£ is supported. However, just 
as for local operators, surface operators can have conformal anomalies which introduce a 
non-trivial scale-dependence. Conformal invariance, however, restricts the functional form 
of the expectation value to be 



where in principle both C and A can depend on the coupling constant of the theoryo. We 
have calculated these quantities in the semiclassical approximation of the gauge theory 
and also in the probe approximation, where a surface operator Oj: is described by a 
configuration of D3 branes. In both regimes we find that 



In the case of local operators one can determine the dimension by calculating a 2-point 
function. We would like to point out, to avoid confusion, that the analogous quantity here is the 
vacuum expectation value of a single sphere, not two. In both cases there is a single length-scale, 
the distance between the points or the radius. 




C = l, 



A = 0. 



(1.2) 
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This result is rather different from the one obtained for surface observables in the six- 
dimensional (0, 2) theory dual to string theory on AdSj x S"^, where the dimension A is 
non-zero [|ig [0 fg] IH [|l3 . 

We then calculate the correlation function between a surface operator Ox; and local 
operators O in Af = 4 SYM. We are able to calculate these correlation functions using 
the three different realizations of surface operators, one using gauge theory, one using 
probe D3-branes and one using the "bubbling" supergravity solutions. We find remarkable 
agreement between these correlators computed in the three different regimes. 

Such two-point functions are obvious quantities to study since local operators serve as 
natural probes of surface operators 0-£, and it is interesting to unravel how they interact 
with them. Moreover, when the surface operator Oj: is supported on an S\ it can be 
approximated at large distances by a linear combination of local operators [TU\. In other 
words, there is an operator product expansion of a surface operator Oj: in terms of local 
ones O. The correlators of the surface operator Oj: with the local operators O determine 
the coefficients in this expansion. 

One particularly interesting local operator is the energy-momentum tensor Tmn- The 
correlator between Tmn and a local operator O is proportional to the conformal dimension 
of the operator O. When we calculate the correlator between Tmn and a surface operator 
Oj: we find a non-zero answer. This is not in contradiction with the already stated result 
that there is no radius a dependence in the vacuum expectation value of a surface operator 
Oj:. As we summarize in in Section 2.2, a surface operator Oj: has three different possible 
conformal anomalies and the two-point function of a surface operator with the energy- 
momentum tensor may calculate a different combination of these three anomaly coefficients. 

We then calculate the correlation function between a surface operator Oj: and Wilson 
and t'Hooft loop operators in A/" = 4 SYM. We are also able to calculate these correlation 
functions using the three different realizations of surface operators, one using gauge theory, 
one using probe D3-branes and one using the "bubbling" supergravity solutions. We 
also find remarkable agreement between these correlators computed in the three different 
regimes. 

The geometrical properties of a correlator between a surface operator Oy, supported 
on a surface E in spacetime and a loop operator supported on a curve C in spacetime 
are quite interesting, since in four dimensions a surface E can be linked by a curve C. 
The correlators of 0-£ with Wilson and 't Hooft loops exhibit rich dependence on all the 
parameters that a surface operator Oj: is characterized by. 
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The Wilson loop we consider was constructed in [15| and preserves eight supercharges. 
This Wilson loop links the surface operator Oj: and depends on two parameters {Oq^i/jq) 
which encode the coupling of the six scalars in A/" = 4 SYM to the loop. It has been argued 
in il^ that the expectation value of this Wilson loop is captured by a matrix model. We 



expect a similar matrix model to be valid also in the case when a surface operator 0-£ is 
inserted into the path integral. This suggests how to extrapolate the value of the correlator 
between O-^ and We^^^^^^ to strong coupling by including the strong coupling answer of the 
matrix model. 

In the realization of the surface operator Oj] in terms of a probe D3-brane, we cal- 
culate the correlator between Oy. and Wg^^^^ by evaluating the on-shell action of a string 
worldsheet that goes between the boundary of AdS^ x and which ends on the probe 
D3-brane in the bulk, and we show that it reproduces the matrix model result. 

In the realization of the surface operator 0-£ in terms of a "bubbling" supergravity so- 
lution, we describe the Wilson loop We^^^^^^ by classical worldsheets in the exact "bubbling" 
supergravity background. Remarkably, the string worldsheet is described by straight line 
in the auxilliary space of the "bubbling" geometry. The correlator between and VFecV'o 
is then obtained by summing over all string worldsheets that are saddle points of the path 
integral, reproducing the matrix model result. 

Analogous calculations are performed for a supersymmetric 't Hooft loop T^q^^q, but 
in this case there is no known matrix model description of the operator and the gauge 
theory analysis is performed in the semiclassical regime. 

The paper is organized according to the different realizations of a surface operator Oj: , 
and in each realization we perform several calculations. In Section 2 we consider the gauge 
theory description of 0-£, where we study conformal anomalies, the vacuum expectation 
value of Ce, as well as the correlators of Ce with local and Wilson and 't Hooft loop 
operators We^^^g and T^q.-^o- 

In Section 3 we study surface operators and their correlation functions using a probe 
D3-brane description. Here we evaluate the classical on-shell action of the D3-brane and 
find the overall scaling anomaly as well as the correlator of Oj: with chiral primary oper- 
ators in = 4 SYM. We also find the string solution describing the Wilson loop We^^^^ 
in the presence of a surface operator Oj: and calculate the the correlator between the two 
operators. We extend this analysis to the case when the Wilson loop is replaced by the 't 
Hooft loop T^o^^Q. 



In Section 4 we turn to the exact supergravity "bubbling" solutions describing surface 
operators. From the exact supergravity solution we find the correlator of O^, with various 
local operators O using holographic renormalization. We then find the exact description 
of the Wilson loop VF^q^^q in terms of a classical worldsheet in the full supergravity back- 
ground corresponding to a surface operator Oj:. We show that these worldsheets have a 
very simple description in the "bubbling" solution and calculate the correlator between Oj: 
and Wgici/iQ by evaluating the on-shell action of the string worldsheet. A similar analysis 
is considered for the 't Hooft loop Tg^^^^ 

We end with a discussion where, in particular, we speculate about a possible matrix 
model dual of a surface operator Oy.- Some computations are relegated to the appendices. 



2. Gauge Theory Description 

In this section we perform calculations with surface operators in A/" = 4 SYM in the 
semiclassical approximation of the gauge theory. We start with a brief review of disorder 
surface operators [|1| (see also 0) and their semiclassical description. We then calculate 
their expectation value, the correlators of a surface operator with local operators and the 
correlators of a surface operator with a Wilson and an 't Hooft loop. 

2.1. Surface Operators in J\f = 4 SYM 

A disorder surface operator 0-£ supported on a surface E C i?^ is characterized by 
a codimension two sing ularity for the AT = 4 SYM fields! In this paper we are inter- 
ested in maximally supersymmetric surface operators in A/" = 4 SYM which preserve the 
maximal subgroupll PSU{1, 1|2) x PSU{1, 1|2) x U{1) C PSU{2, 2|4) of the four dimen- 
sional superconformal group of A/" = 4 SYM. The corresponding Euclidean supergroup 
has an 5'0(1,3) x SO{2) x SO{4) bosonic subgroup, of which the first component acts on 
space-time alone, the last is an /^-symmetry and the 5*0 (2) is a diagonal combination of 

^ For previous work involving codimension two singularities in gauge theory see [p^ [|17|. 
^ This is the supergroup for a maximally supersymmetric surface operator in R^'^. In this 
paper we consider the theory in i?^, but surface operators exist in both. 
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space-time and i?-symmetries. S'0(l,3) is the Euclidean conformal group in two dimen- 
sions and the only 2-dimensional submanifolds of which posses this symmetry are 
and S'^, which will be the locus of support of our surface operators!. 

The symmetries allow a complex scalar field in the Af = 4 SYM multiplet to develop 
a pole near E, reducing the S'0(l,5) x SU{4) symmetry of A/" = 4 SYM to 5'0(1,3) x 
SO{4) X [/(I). 

An alternative way to study a surface operator is to consider the gauge theory oni 
AdSs X 5"^ instead of R^, the two geometries being related by a Weyl transformation. Then 
a surface operator 0-£ is characterized by a non-trivial boundary condition for the fields 
in the gauge theory near the boundary of AdSs x . The advantage of this formulation is 
that the symmetries of the surface operator 0-£ are made manifest, as they are realized as 
isometries of AdS^ x instead of as conformal symmetries in R^. The choice of surface 
E is now captured by the topology of the AdS^ x 5"^ conformal boundary. For T, = S'^ one 
must consider the AdSs metric in global coordinates while for E = i?^ one must consider 
the metric in Poincare coordinates. 

A surface operator in i?^ induces a non-abelian vortex configuration for the 
gauge field near the surface operator C^;, which produces a codimension two singular- 
ity at E. The gauge field configuration breaks the gauge group G to a subgroup L C G 
near E. The singularity in the gauge field corresponding to a surface operator Ox; with 
L = Ufii U{Ni) C U{N) is given by 



A 



/ctiOlAT, ■■■ \ 
a2 0liV2 ■■■ 

V ■■■ aM ® Inm ' 



dtl:, (2.1) 



where ij) is the polar angle in the local normal bundle to E and 1^ is the n-dimensional 
unit matrix. The parameters ai take values in the maximal torus — /Z^ of the 
gauge group G = U{N). 



^ The two surfaces are related to each other by a global conformal transformation. On R^'^ 
there are also hyperboloids with geometry H2 and dS2 which are conformal to space-like and 
time-like planes respectively. One should be able to construct maximally supersymmetric surface 
operators on them too. 

^ Throughout this paper AdS refers always to Euclidean AdS space. 
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The operator Oj: also introduces into the path integral two-dimensional ^-angles that 
couple to the unbroken t/(l)'s along S. This can be represented by the operator insertion 



exp 




(2.2) 



where F^''^ is the field-strength for each of the unbroken U{1) factors. The parameters 
rji take values in the maximal torus of the 5'-dual or Langlands dual gauge group 
Since the Langlands dual of G = U{N) is also = U{N)^ we have that the matrix of 
^-angles of a surface operator Oj: is also given by an L-invariant matrix 



^771 ^Iat, 

??2 ® liVa 








(2.3) 



A maximally supersymmetric surface operator 0-£ also excites a complex scalar field 



in A/" = 4 SYM which we label by $ 
E of the form 



71^ 



-f It develops an L-invariant pole near 



/32 -I- i72 ® 1 N2 



V2; 








(2.4) 



where z 



re 



lip 



is the complex coordinate in the local normal bundle to E. This form 



of the singularity is completely determined by imposing superconformal invariance and 
L-invariance. 

In summary, a maximally supersymmetric surface operator O-^ in A/" = 4 SYM with 
Levi group L = YltLi U{Ni) is labeled by AM L-invariant parameters (ct^, 7^, 77^). The 
surface operator in i?^ is defined by performing the Af = 4 SYM path integral over smooth 
field fiuctuations around the L-invariant singularities (|2.1| )( p7^ ) with the insertion of the 
operator ( p.2|) . We must mod out the path integral by gauge transformations that take 
values in L C U{N) when restricted to E. The surface operator 0-£ becomes singular 
whenever the parameters that label the surface operator {ai, f^i^ji^rji) for / = 1,---,M 
are such that they are invariant under a larger symmetry than L, the group of gauge 
transformations we have to mode out by when evaluating the path integral. 
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In the AdSs x picture, the surface operator Oj: is described by a non-trivial holon- 
omy of the gauge field ( |2.1|) around the non-contractible 5"^ parametrized by the coordinate 
"0 and by the insertion of ( |2.2[ ), where E is now the conformal boundary of AdSs x S^. 
The operator 0-£ is also characterized by a non-singular configuration for the scalar field 



71" 



f Pi + hi® In, •■■ ^^ 



V ■■■ Pm + ^m^Inm' 



(2.5) 



where is the coordinate parametrizing the in AdSs x S^. Therefore, a surface operator 
is described by a half BPS vacuum state of the gauge theory on AdSs x S^. 
Note that the field configurations for E = i?^ and for S'^ are exactly the same when the 
gauge theory is on AdSs ^ S^. For A/" = 4 SYM on on the other hand, the description 
of the singularities of the fields produced by a surface operator with E = S"^ of radius 
a is quite different from E = i?^. One can perform a Weyl transformation from global 
AdSs X S^ to i?* and use that the gauge field and the scalar fields transform as weight zero 
and weight one fields respectively. Starting with the metric of in spherical coordinates 

ds^ = dr^ + r^id'd'^ + sin^ 'd dip^) + dx^ , (2.6) 

and defining 

r = —5 = — TT , r = r smh p , x = r sm ip , 

4a^ (cosh p — cos 'i/')^ 

(2.7) 

we find the metric 

ds^ = + ginh^ p (d^^ + sin^ ^ dif^) + rfV^) , (2.8) 

which is conformal to AdSs ^ iii global coordinates. 

In the coordinate system ( |2.6| ), a surface operator Oy, with E = S"^ C -R^ located at 
r = a produces a singularity on the scalar fields given by ( p.4|) with the replacement of 
z — > re*"^, where i/j is written in terms of the coordinates ( |2.6| ) using ( |2.7| ). The gauge 
field singularity is exactly the same as before ( |2.1|) , where again di(j is written in terms 
of the coordinates ( p.6| ) using (|2.7|) . Finally, the insertion of the ^-angles ( p.2|) is now 



obtained by integrating TrF; over the 5"^ parametrized by ^? and f in (p.6|) at r = a. 



9 



2.2. Surface Operator Expectation Value 

We now compute the expectation value of the surface operator O^, in the semiclassical 
approximation of the gauge theory. We want to calculate the expectation value for both 
E = i?-^ and T, = S"^. When the operator is supported on we anticipate that the 
expectation value will be one by supersymmetry. For the case T, = S"^ the expectation 
value may depend non-trivially on the radius a of the S'^. 

Since we are dealing with operators defined on a surface E, an even dimensional sub- 
manifold in spacetime, the surface operator Oy: may posses conformal anomalies [p!0| [p!T|. 
The conformal anomalies for surface operators are given by local expressions that depend 
on the puUback onto the surface E of the curvature invariants in spacetime and also on 
the extrinsic curvature, which encodes the embedding of the surface E in spacetime. 

The classification of the surface conformal anomalies can be rephrased as a cohomology 
problem, as in the case of usual conformal anomalies [^. One must classify solutions to 
the Wess-Zumino consistency conditions modulo terms that can be obtained as the Weyl 
variation of a local term. A basis of the Weyl cohomology has been recently presented 
in (for related previous work on surface anomalies in six dimensions see [p!2| ||13|| [[I^ ). 
The surface anomaly depends on three anomaly coefficient c^, i = 1, 2, 3 

Aj: = ^J^ d'aVh (^cii?(2) + C29mn [h^"h^'' - \h^^hP"^ K^.K^, + c^.h^'' h^^ C , 

(2.9) 

where h^^, and Qmn denote the induced metric on E and the spacetime metric respectively 
while Cmnpq and denote the Weyl tensor and the extrinsic curvature. The anomaly 
coefficients Ci are the analogs of the well known a and c anomaly coefficients in four 
dimensions. The expression for ci depend on the choice of field theory as well as on the 
choice of surface operator in that field theory. 

For the case E = i?^ C -R^ we have that Aji2 = 0. For E = S"^ C -R^ one gets 

^ ^ V^i?(2) ^ 2 , g^^ (^^^'^h^P - h'^'-hP'^^ KJZK"^, = Cpp^a = 0. (2.10) 

Therefore As2 = 2ci. 

Under a conformal transformation Qmn ^'^'^dmn the expectation value of a surface 
operator Oj: supported on S"^ C changes by an amount proportional to the anomaly 

5,(Oe) = 2ci(0e). (2.11) 
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Therefore the expectation value of 0-£ is determined in terms of the radius a of the 5"^ by 



(C»E)oca^^i. (2.12) 

The precise form of the expectation value of Oj: now depends on the chosen operator and 
on the theory under discussion, which uniquely determine the anomaly coefficients c^. 

We now focus on our specific case of the maximally supersymmetric surface opera- 
tors Oj: in a/" = 4 SYM and calculate their vacuum expectation value in the semiclassical 
approximation. This is achieved by evaluating the classical SYM action on the field con- 
figuration produced by the operator 

(Ce) =exp(-5)|surface. (2.13) 

It is easiest to perform this calculation by considering the field configuration produced by 
a surface operator Oj: on AdSs x S^. The relevant part of the action of A/" = 4 SYM on 
AdSs X is given by 

JZ= (\D^\^ + ^\^ A , (2.14) 



9ym " V' " ' ' 6 

where is the conformal coupling for the scalars and R^'^^ = —6 is the scalar 

curvature of AdSs x 5"^. The surface operator Oj: induces a non-trivial field configuration 
on the field $ (|2.5|) on AdSs x which depends only on the angular coordinate ifj of the 
through $ = $0^"*^, where $o is a constant diagonal matrix ( |2.5| ). The operator Oy. 



also induces a nontrivial gauge field configuration (|2.1| ) which nevertheless commutes with 
$0? so that = 9$. Therefore = and the two terms in the Lagrangian ( |2.14| ) 

exactly cancel each, so that 

(Oe) = 1. (2.15) 

This result holds regardless of whether we use the Poincare patch or global AdSs metric 
in the field theory, so that (Oe) = 1 for both E = i?^ and T, = S'^. 

The same result can also be derived by studying Af = 4 SYM on R^. In this case the 
scalars do not have a mass term. However, one must include in the action a boundary term 
along the surface operator Oe that guaranties that the boundary terms in the equations 
of motion vanish. To see that, consider the relevant part of the action in fiat space 

S=^ [ d^x^g"^^ Tr [0^^ D^^) (2.16) 
9ym J 
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with the metric (|2.8|) . In deriving the equations of motion there is a boundary term near 
the locus of the surface operator at a cutoff p = po 

5S= ^— [ (I'd d<p dtp f"^ sinh^ pQsin'd {5^ dp^ + 5^ dp^) . (2.17) 
9ym J 

This does not vanish, since $ diverges. To fix this we need to add to the action the 
boundary term 

■^boundary = 2 / ^'^ ~ ^^^^^ PO ^ ^r ($$) . (2.18) 

9ym J ^ 

The variation of the bulk and boundary action would now give the boundary equation 
9p$ — r sinhpo $/a = 0, which is indeed satisfied by our solution. 

Evaluating the classical action with a cutoff po gives the divergent result 

S=^l dpdtPd^d^ ^"^^f ^"^'^^^"^ f;iVKA^+7f) = 4vr^inh^pof:iVKA^+7f) ■ 
9ym J coshp-cosV' f-[ 

(2.19) 

The boundary term ( |2.18|) , though, exactly cancels this bulk piece leading to a total action 
that vanishes, exactly as we found in the AdS'i x 5"^ picture. 

It follows from ( |2.12| ) and from (Oe) = 1 for a spherical surface operator Oy. in A/" = 4 
SYM that the anomaly coefficient ci for this theory (and this family of surface operators) 
is ci = 0. It would be interesting to calculate separately the anomaly coefficients C2,C3, 
which can be done by computing the anomaly for more general surfaces E as well as by 
considering A/" = 4 SYM in more general curved four dimensional backgrounds. 

2. 3. Correlator of Surface Operators with Local Operators 

We now proceed to calculate the correlation function of a surface operator Oy, with 
various local operators in A/" = 4 SYM in R^. We consider the correlator with the dimension 
A = 2, 3 chiral primary operators (CPO's) 0/\ of A/" = 4 SYM, the traceless stress-energy 
tensor Tmn and with the U{1) spacetime symmetry current which generates rotations 
in the space transverse to the surface E. 

In the semiclassical approximation, the correlation function of a surface operator Oy, 
with a local operator O is obtained by evaluating O in the background field produced by 
the surface operator Oy 

{Oy ■ O) 



{Oy) 
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C^lsurface- (2.20) 



If instead we consider the surface operator Oj^ as a non-trivial boundary condition for 
A/" = 4 SYM on AdSs x , the correlation function of 0-£ with a local operator O in 
corresponds to the vacuum expectation value of the local operator O in the state |Ce)- 
Thus on AdSs x S'^ the correlation function ( pj.20| ) in R"^ corresponds to 

(^Oo.- (2-21) 

In all these computations the position dependence of the local operator is determined by 
conformal Ward-Takahashi identities. In the case of E = i?^ C the correlator of a 
surface operator 0-£ with a local operator O of dimension A scales with the distance r 
between the local operator O and the plane where the surface operator Oj: is supported 
byB 

• Co (2 22) 

(Oe) "r2A- ^2.22J 

The quantity that we need to compute for the various local operators O is Cq. 

When the surface operator O-e is supported on 5"^ C i?^, the correlator is given by 

{O ■ Oe) Co 



where f is the conformally invariant distance (|2.7|) 



(2.23) 



X 



2 



(2.24) 



2a 

combining a, the radius of the 5^, the radial position of the local operator in the 
including the surface, r and the transverse distance to the local operator, x. In principle 
the proportionality constant in (p.22| ) (|2.23|) Co should be the same for = S"^ and 
E = i?^, except that there is the possibility that some contributions that exist for the 
sphere are "pushed to infinity" in the case of the plane and will not show up. This does 
not happen in the specific calculations we performed but may show up once one includes 
quantum corrections. 

When one considers the gauge theory on AdS^ x 5"^ the surface operator Oy. is sup- 
ported at the boundary of AdS^. Therefore, there is no invariant distance like r or f in 
i?^. In this case 

{0)o^=Co. (2.25) 



^ For operators that are not scalars, the tensor structure of the correlator is determined by 
Lorentz invariance. 
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The CPO's are scalar operators transforming in the (0, A, 0) representation of the 
SU{4) R-symmetry. The operators are given by 



Oi = '^^0^Ci„,^Tr . . . f-) , (2.26) 



where A = Qym^ Hooft couphng and = C-_^ ^^x^^ ...x'^^ are the SO{6) 

scalar spherical harmonics. The operators ( |2.26p are normalized such that their two point 
function is unit normalized 

{Oi{x)Oi{y)) = (2.27) 
\x y\ 

A surface operator 0-£ preserves an SO{4) subgroup of the SU(4) R-symmetry group. 
Therefore, the correlation function of a chiral primary operator (9^ with a surface operator 
0-£ is nontrivial only if the CPO is an SO{4) singlet. By looking at the decomposition of 
the (0, A, 0) representation of 5*0(6) ~ SU{4) under the obvious SO{4) x 50(2) subgroup 
one finds that there are A + 1 CPO's of dimension A which are singlets under 50(4). 
These 50(4) invariant operators have charges k = —A, — A + 2, • ■ ■ , A — 2, A under 50(2). 
We label the 50(4) invariant spherical harmonics by Y^'^ = O^^'^j^x*^ . . .x^'^ (for more 
details see Appendix A). With these spherical harmonics we can write down the 50(4) 
invariant CPO's as 

Oa,, = J^^CtXT^ (2.28) 

Using the explicit expression for the 50(4) invariant spherical harmonics with A = 
2, 3, we can write down the unit normalized CPO's 



Os, = ^Tr (2*^4 - ; O3.3 = 



(2.29) 



where as before $ = -^(^ + i(p ) and 0a, -fc = 0A,fc- 
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The semiclassical correlations function ( |2.2C1| ) of these CPO's with the surface operator 
0-£ are then given by 

(02.0 ■ 0e) 1 StT^ ^ 2 2 (02,2 ■ ^e) ^ ^ m ( R ^ ' ^2 



(03,3 -Oe) ^ - - 



(2.30) 

The remaining correlators can be obtained by complex conjugating ( |2.3C1| ). 

These result are purely classical and they may receive quantum corrections. The 
surface operator Oz breaks near E the gauge group from U{N) to L = UZi U{Ni), and 
also breaks /^-symmetry group down to SO [4). This breaking may introduce different 
interactions among the different fields at the quantum level. We therefore expect the 
quantum corrections to depend on the choice of Levi group given by the rank of the 
unbroken gauge groups A''!, ■ ■ ■ , Nm- The calculations performed using probe D3-branes 
and the "bubbling" supergravity solution in the next two sections indicate, however, that 
the loop corrections terminate after a finite number of loops! The perturbative corrections 
to the correlator of O^, with OA,k terminate at ordei^ ^(^-l'=l)/2 n would be interesting 
to understand the truncation of the loop corrections directly in the gauge theory. 

We now consider the correlator of Oj: with the stress-energy tensor of A/" = 4 
SYM in R^, which was already computed in 0. In the presence of a planar surface operator 
in the correlator is fixed up to function h 

^^^^ = ft^. = (r„<-O.> = 0. (2.31) 

Here x"^ = (a;^,a:*), where x'^ are coordinates along T, = R^ and n* = x'^/r is the unit 
normal vector to the planar surface. If A/" = 4 SYM is on AdSs x the expectation value 
of the stress-energy tensor Tmn in the presence of a surface operator Oj: is given by 

{Tab)oj: ^ h 9ab {T^^ip) q = " 3/i, (2.32) 



^ Note that A — |fc| is always even, so we get a standard perturbation expansion in the 't Hooft 
coupling A. 
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where Qab is the metric on AdSs and i/j is the coordinate on the 5"^. 

The piece of the stress-energy tensor of A/" = 4 SYM which involves the gauge field and 
the complex scalar field $ excited by the surface operator Oj: is given by (see e.g. pPpi 



Tr 



+ 



9ym 



Tr 



(2.33) 



Using the semiclassical formula ( |2.2C1| ) we find that 



h 



M 



1=1 



(2.34) 



or equivalent ly 



M 



(Tat) 



'^9ym 

M 



(2.35) 



9ym 



7z 



1=1 



In this case, we expect that the only perturbative correction appears at the one loop level, 
of order A, which as we will see is captured by the "bubbling" supergravity solution. 



In the analogous coefficient h for Wilson and 't Hooft loop operators was dubbed 
the scaling weight, the name suggesting that it should generalize the notion of conformal 
dimension of local conformal fields to extended objects. The analogy with conformal 
dimension of a local operators is not complete. The conformal dimension of a local operator 
measures its behavior under conformal transformations, while for surface operators there 
are three different anomalies under such transformations which we dubbed Ci in ( |2.9| ). We 
expect the scaling weight h to be related to the anomaly coefficients Ci and it would be 
interesting to find the explicit relation. 

We also calculate the correlator of Oj: with the conserved current J^. This current 
generates the U{1) symmetry that acts by shifts on the 5"^ of the AdS^ x geometry on 
which the gauge theory is defined. However, it is not a symmetry of the surface operator 
as the singularity induced on the scalar field (|2.4|) is not U{1) invariant. In conformal field 
theories, the current associated with a spacetime symmetry is constructed by contracting 



We neglect the scalar potential term, as it does not contribute to the correlator. 
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the stress-energy tensor Tmn with the conformal KiUing vector generating the symmetry 
im = Tmnia- Therefore, the correlator of with O-^ is given by 



Jt),, ={T^^)o. = ;^J2Ni{Pf + ^f). (2.36) 
' C's 9ym 

2.4- Correlator with Wilson and 't Hooft Loop Operator 

The last calculation we perform is the correlator of a surface operator 0-£ with a 
Wilson loop operator and also make comments about the correlator of a surface operator 
Oy: with an 't Hooft loop operator. This combination of operators is particularly interesting 
since in four dimensional space a two dimensional surface can be linked by a curve. 

We take the Wilson loop supported on a curve C to link the surface E on which the 
surface operator Oj: is supported. The simplest way to realize this configuration is by 
taking the surface T, = C and consider a circular loop C = S"^ in the orthogonal 
plane. A natural choice of Wilson loop to consider is the 1/4 BPS circular Wilson loop 
with periodic coupling to the scalar fields with arbitrary phase shift i/^o [H 



W^^ = j^Tr exp J {iA^ + r cos(i/' - tPo) (j)^ - r sin(i/; - ^o) cj)^) dtP . (2.37) 



Here (j) and (j) are the real and imaginary part of the complex field turned on by the 
surface operator (9^, given by $ = -^{'P^ + i4>^)- 



We consider an even more general Wilson loop operator, presented in [15| and studied 



further in [^, with an extra coupling to a third real scalar field (f)^ in the J\f — 4 SYM 



multiplet, with a relative coupling given by the angle 9o 

6*0 (/>^ + ^2 sin6'o3f?(;2$e-^'^o)) . (2.38) 

Fort 6*0 = 7r/2 it reduces to ( |2.37| ) (note that the ip dependence comes now from z = re*"^), 
while for Oq = it's the well studied maximally supersymmetric circular loop ||24]| [p5 |. 

Note that this Wilson loop operator depends on two parameters 9o and ipo and op- 
erators with different values of these parameters are quite different. In particular, the 
supercharges that are preserved by the Wilson loop Wq^^^^ depend on these parameters, 
though they are always compatible with those of the surface operator. This is proven in 
Appendix B. 
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Evaluating the correlator of a surface operator Oj: with the Wilson loop (|2.38| ) in the 
semiclassical approximation amounts to evaluating the Wilson loop operator ( pj.38| ) on the 
fields produced by the surface operator (|2.1J ) (p^). The result is given by 

M 



class 



1=1 



27r sin 9o{Pi cos tpo + 'ji sin ipo) + 2niai) 



(2.39) 



The extra scalar cj)^ does not contribute at the classical level, as the surface operator O-e 
does not excite it. 



In the absence of a surface operator O^, it is conjectured in that the exact 
expectation value of the Wilson loop Wg^^ip^ is given by summing up all ladder diagrams 
in the Feynman gauge, which is captured by a matrix model. In the limit of large 't Hooft 



coupling A it has an asymptotic expansion with two saddle points which can be written 
schematically as 



{Weo,i^o) ^exp 



COS 60 



+ exp 



cos 9o 



(2.40) 



This result, including both saddle points can be reproduced from string theory. 

The expression ( |2.4C1| ) can be regarded as coming purely from the scalar field (p'^ , which 
appears with a factor of cos^o in the loop ( p.38|) , while the other two scalars are present 
to guarantee supersymmetry. 

The presence of the surface operator can, of course, modify the contribution of (f)^ to 



the expectation value of the Wilson loop ( |2.38| ) , but we now argue that its contribution is 
well under control. Indeed this is what happens for the Wilson loop VFso-V'o presence 
of a chiral primary local operator (see also [^[|2^), where the exact result is given 
by a normal matrix model, but the asymptotics are still governed by an exponent like in 

In particular, one possible effect of the surface operator O-^ on the scalar (p^ in the Wil- 
son loop (p.38| ) is in breaking of the U (N) gauge group to the Levi group L = nf=i ^ (-^Oi 
thus reducing the number of degrees of freedom in the system. Now, in each U{Ni) sector, 
(p^ gives the same matrix model result with the replacement A —>■ XN^/N. Assuming that 
all the Ni are still large, so that we may still use the approximation in (|2.40|) , we predict 
that the correlator of the Wilson loop W^^^-tp^ with a surface operator Oy, at strong coupling 
is given by 

M 




exp 



COS 9q 



+ exp 



cos 6*0 



(2.41) 



X exp [27rsin^o(A cosi/'o + 7z sinipo) + 27Tiai)] . 
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We will later see that this result is reproduced by performing a string theory calculation 
in the "bubbling" supergravity background dual to a surface operator. 

We note that unlike local operators, this class of Wilson loops are sensitive to the 
Aharonov-Bohm phase ( |2.1|) induced by the surface operator Oj:. 

Similarly to the Wilson loop, one can consider an 't Hooft loop, which should yield 
the S-dual result (accounting of course also for the S-duality transformation of the surface 
operator parameters («/, 7;, 77;)). Unlike the Wilson loop, the 't Hooft loop does not 
have a simple weak-coupling description. Rather, much like the surface operator Oj:, it 
can be thought of as requiring certain singularities for the gauge field and scalar field along 
the loop. 

The 't Hooft loop T^q^^q we want to consider should be the dual of the Wilson loop 
( p.38|) , which couples to three scalars and depends on the parameters 9o and i/'o- An 't 
Hooft loop can be represented by a Dirac monopole embedded along some U{1) subgroup 
of U{N) with the accompanying scalar source. The field-strength around the monopole 
has the form 



where 9 and (j) are local coordinates on a sphere surrounding the loop, and is a generator 
in the Cartan subalgebra of U{N). This should be the behavior of the gauge field close 
to the 't Hooft loop Tg^^^^^ for the straight line. In the classical approximation this is the 
exact solution. So consider now the spherical surface operator Oj: wrapping the straight 
't Hooft loop Tsq^^q. The surface operator couples to the magnetic field produced by the 
't Hooft operator through the parameters 77 in (|2.2|) entering the definition of (9e 



This produces the S'-dual of the phase exp(27rzQ;) in (|2.41|) obtained from the correlator of 
Oj: with the Wilson loop ( |2.38| ), as a ^ 77 under ^'-duality 

The contribution from the scalars $ and $ is almost identical to the one in the Wilson 
loop calculation, as the effect of the singularity on the scalar fields produced by the 't Hooft 
loop can be modeled in the semiclassical approximation by inserting into the A/" = 4 SYM 
path integral the operator 



2 



(2.42) 




(2.43) 




(2.44) 
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Inserting the field ( |2.4|) produced by the surface operator Oy: in (|2.44| ) gives the 5'-dual 
of the Wilson loop result ( ^T^ ) as \P + i'fl -> 4n/g^j^\(3 + i'y\ under ^'-duality The 
contribution of 0^ as in (|2.4C1| ) has so far not been reproduced from the 't Hooft operator in 
a gauge-theory calculation even in the absence of the surface operator Oj:, and we do not 
attempt to recover that here. We will see in Section 3 and in Section 4, where we calculate 
the Wilson loop and 't Hooft loop using string theory, that the resulting expressions are 
indeed related to each-other by S'-duality. 

We note that the correlator of Oj: with an 't Hooft operator is sensitive to the two 



dimensional ^-angles (|2.2|) that appear in the definition of a surface operator O^. 



3. Probe Description in AdS^ x 

In this section we perform calculations with surface operators in A/" = 4 SYM in the 
probe approximation of the dual Type IIB string theory on AdS^ x S^. We start with a 
description of surface operators in A/" = 4 SYM in terms of D3-branes in AdS^ x S^. We 
then calculate their expectation value, the correlators of a surface operator with CPO's 
and the correlators of a surface operator with a Wilson and an 't Hooft loop. 

3.1. Surface Operator as D-branes in AdS^ x 5"^ 

A surface operator Oj: in A/" = 4 SYM corresponds in the probe picture to having 
a D3-brane in AdS^ x ending on the boundary along the prescribed surface E. This 
D3-brane solution in AdS^ x was first studied by [p9| . 

We start by describing the details of the AdSs x D3-brane solution that we need 
to perform the various calculations in this section. We take the metric on the to be 

dn5 ^ cos'^ edfls + de^ + sm^ed(l)'^ , (3.1) 

which makes manifest the SO{4) x SO{2) C 5*0(6) symmetry of the S^. For the surface 
operator Oy, supported on Tj — it is convenient to write the bulk AdS metric in two 
different slicings 

dsf^-^ = ^ {dy'^ + dxl -t- dx\ + dr"^ -f dip"^) , (3.2) 

(is(2) = "^^^ ^ ^ {dz^ -f dx^ -f dx\) -\- du^ -f sinh^ udtp'^^ (3-3) 
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where we work in units where the AdS^ and radius of curvature is one so that = 
AnQsNl^ = 1. The two coordinate systems are related by 

2: 

y= — - — , r = ztanhw. (3-4) 
cosh u 

The first choice of coordinates ( p.2|) is more suitable when considering the dual A/" = 4 SYM 
theory on i?^. The second one (|3.3|) is more suitable when A/" = 4 SYM is on AdS^ x S^, 
where AdS2, x is the metric on the conformal boundary located at w ^ 00. 

We choose the RR four-form potential in the analog of the Fefferman-Graham gauge 
||30|| for the metric, so that the four-form has no components transverse to the AdS^ 



boundary. Therefore, we take the RR four-form potential in the metric (|3.2|) and (|3.3|) to 
be given by 



C^^^ = dxo A dxi A dr A di(^ , (3-5) 



r 



y 



C\ ' = T — dz A dxo A dxi A dtp . (3-6) 

Z"^ 

The worldvolume coordinates of the D3-brane embedding in the coordinate system 
( p.2|) are given by xq, xi, r and tp, and the D3-brane has non-trivial embedding functions 
y = y{r) and on the 5"^ in the coordinate system (|3.1|): 6 = ^/2 and (p = </>('(/')• With this 
ansatz the D3-brane action is given by 

Sd^ = Td3 j dxQ dxi dr dip ^ + + y^<p'^) - T^s j dxQ dxi dr dip ^ (3.7) 

where the D3-brane tension is given by Tjy^ = (in the units where = 1). The 
half-BPS solution is given by |^ 

y{r) = -r, ip + (p = (po , (3.8) 

K 

where and (po are integration constants. 

In the coordinate system ( |3.3| ) the D3-brane solution has worldvolume coordinates z, 
x*^, x^ and tp and the embedding is given by 6* = 7r/2 in ( |3J] ) and 

sinh uq = K , ip + (p = (po . (3-9) 

The induced metric on the brane is that of AdSs x 5"^, so that we have the freedom of 
turning on a Wilson line for the gauge field A and dual gauge field A living on the D3-brane 
worldvolume along the non-contractible S^. Therefore, the D3- brane solution ending on 



cosh"* u 
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E = i?^ on the boundary depends on four parameters which can be identified with 
those of the gauge theory description of a surface operator O^, via0 

(3 + ij = ^smhuoe'^% (3.10) 

Therefore, in the probe approximation a surface operator 0-£ with unbroken gauge 
group L = YliLi U{Nl) and labeled by the 4M parameters {ai, /3i,'ji,rii) corresponds to 
a configuration of M stacks of D3 branes in AdS^ x 5"^ ending on the boundary on the 
surface E, with each stack consisting of A^^ coincident D3 branes. In Appendix C we study 
the supersymmetry preserved by the probe D3-brane brane and show that it is the same 
as the supersymmetry preserved by the surface operator 0-£ in A/" = 4 SYM. 

For the probe D3-brane description of the surface operator Oj: in A/" = 4 SYM sup- 
ported on E = S"^ it is convenient to write the bulk AdS^ metric in two different slicings 

dsf^^-j = ^ {dy"^ + dr^ + r^{d^^ + sm'^^d^^) +dx^) , (3.11) 

dsf^) = cosh^ u(rfp^ + sinh^ p {d-d"^ + sin^ 'ddip'^)) + dv^ + sinh^ u dif)"^ . (3.12) 
The two coordinate systems are related by 

y = — -, r = J/ cosh u sinh p , a; = y sinhw sin?/', (3.13) 

cosh u cosh p — sinh u cos ip 

where a is the radius of the E = S"^ on which the dual surface operator Oy, is supported. 
Again, the choice of coordinates ( p.ll|) is appropriate when considering the gauge theory 
on R'^ while the one in ( |3.12D is appropriate when the gauge theory is on AdS^, x S^. 
The associated expressions for the RR four-form potential are given by: 

cf^ = — —dr^d{}^dip^dx, (3.14) 
y 

cl^"* = cosh'* u sinh^ p sin dp A d'd A d(p A difj . (3.15) 



Note that in the units = 1, we can write alternatively P + ij = sinh no e"^° (see also 
equation (4.9). 
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The D3-brane worldvolume coordinates in the metric (|3.12|) are p, 99 and tl^, and 
the embedding is again given by 6* = 7r/2 in (|3.1|) and: 



sinhwo = K, i/j + (p = (pQ . (3.16) 

In the coordinate system ( p.ll| ) the D3-brane embedding translates into 

(r2 + x2+2/2_a2)2^4aV ^^2 ^3^7^ 



The map between the parameters of the D3-brane solution and the parameters defining the 
dual gauge theory surface operator Oj: is also given by ( |3.1U| ). As explained earlier, con- 
sidering M stacks of D3 branes yields a surface operator Oy, with parameters {ai, (3i,'yi,r]i) 
for / = 1,...,M. 

3.2. Expectation Value in the Probe Approximation 

The expectation value of the surface operator 0-£ in A/" = 4 SYM corresponding to a 
single probe D3-brane is given in the leading >> 1 and A >> 1 approximation by 

This corresponds to a surface operator 0-£ with unbroken gauge group L = U{1) x SU{N — 
1) and with parameters (a, /3, 7,77). Therefore, we have to evaluate the classical action of 
the D3-brane solutions in the previous subsection. 

For the case of a surface operator O-^ on E = i?^, we find that the action of the dual 
D3-brane in the metric ( p.2| ) is given by 



N r 1 
<Sd3 = TTT / dxi dr dip 



(1 + 2/'2)(^2 _^ ^2^2) _ J. 



(3.19) 



where Sd3 = Sdbi — <Swz- In order to evaluate the classical action it is not enough to 
plug in the solution (/> = ~1 and y = r/n into the action. Examining the variation of the 
action one finds a boundary term that needs to be cancelled. That is due to the fact that 
y' is a non-zero constant and to get a good variational problem it is necessary to subtract 
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Therefore, the total D3-brane action is 



<Sd3 



N 
2^ 



dxo dxi dr dip — 



2A2 



1 + y 



12 



(3.21) 



Plugging in the D3-brane solution (|3.8| ) one sees that the integrand vanishes iSoslon-sheii = 
0. 

It is illustrative to repeat the calculation of the D3-brane action in the second coor- 
dinate system ( |3.3| ), where now the field u = uq is a. constant, so the variational principle 
for it is well defined. The classical on-shell action for the D3-brane is now given by 



S 



D3 \ on — shell 



N 
2^ 



dxo dxi dz di/j 



cosh uq 



cosh Uq \ sinh uq^- (fP — cosh uq 



= 0. 
(3.22) 

We conclude that the expectation value of a surface operator supported on E = i?^ in 
the probe approximation is given by 



(3.23) 



which agrees with the semi classical gauge theory computation (|2.15| ). 

For the case of a surface operator Oy. supported on E = S"^, we find that the action 
of the dual D3-brane in the metric ( p.l2|) is given by 



S 



D3 



N 
2^ 



dp dd d(p dip cosh uq sinh p sin 



sinh^ Uq + (p"^ — cosh uq 



(3.24) 



In this case u = uq is also a constant and it satisfies the proper variational principle 
so after setting ^ = — 1 we find again that the on-shell action also vanishes Sds = 0. 
Thus the expectation value of a surface operator Oj: supported on E = S"^ in the probe 
approximation is given by 

(Oe) = 1, (3.25) 

agreeing with the semiclassical gauge theory computation ( |2.15| ). 

The conformal anomaly (|2.9|) for the surface operator with T, = and S"^ vanishes 
in the probe approximation and due to ( |2.12| ) we have that ci = 0. It would be interesting 
to consider more general D3-brane solutions that end on a boundary with a different metric 
and on a general surface E. This would allows to compute the other anomaly coefficients 
c^ in (U). 
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3.3. Correlator with Local Operators in Probe Approximation 

We now compute in the probe approximation the correlation function of a surface 
operator and the CPO's Oj^ of dimension A = 2, 3 in A/" = 4 SYM. Analogous com- 



putations have been done for Wilson loops in [|T^[^T|[^ and for Wilson surfaces in six 
dimensions in \W. 



A CPO Oj^ in A/" = 4 SYM couples to the five dimensional scalar fluctuation 



3J]. This fluctuation mode diagonalizes the linearized equations of motion of Type IIB 
supergravity expanded around AdS^ x S^. It solves the Klein-Gordon equation in AdS^ 

m 

V^V^/ = A(A -4)s^ (3.26) 

where A is the dimension of the dual CPO (9^. 

In the gauge theory (Section 2.3) we saw that the only chiral primaries that have 
non-trivial coupling to the surface operators are the SO{4) singlets. This can be seen here 
from the fact that in the metric (|3.1|) the brane is at ^ = 7r/2, where all the other spherical 



harmonics vanish. We can combine all the which couple to the probe brane into one 
field on AdS^ x as 

s = Y: E Y^'>^{e=l,<p)s^'\ (3.27) 

A fc=-A,-A+2,--- 

where Y^''^ are the 50(4) invariant spherical harmonics of SO{6) (see Appendix A). When 
evaluated on the D3-brane they are given by 

y^''^(^ = |,</)) =CA,ke'''t (3.28) 

The correlator of a surface operator Oj: with is captured in the probe approxi- 
mation by the solution of the linearized equations of motion of Type IIB supergravity for 
in the presence of the D3-brane source. In order to calculate this we first need the bulk 
to boundary propagator in AdS^. Let us consider at a position (xq, x'l, r', i/;') on the 
boundary of AdS^ in the coordinate system ( |3.2| ). The propagator from the insertion point 
to a point in the bulk is given by 

G — c - , (3.29) 

{{x — x')^ + + — 2rr' cos('i/' — V'') + v'^)^ ' 

25 



where c = ^2-A/2"^y^ chosen such that the bulk computation of the two-point 
function of the dual CPO operator (9^ is unit normalized as in ( |2.27|) . We find it convenient 



to work in the coordinate system ( |3.3| ) where AdS^ is foliated by AdSs x slices. If we 
place at {xq = x'l = 0, ifj' , z' = d) then the bulk to boundary propagator reads 

G = c ^ . (3.30) 

cosh uo{xq + xf + z"^ — 2dzta,nhuo cos{ip — + d"^)^ 

While the bulk field has a simple propagator, it has rather complicated couplings 
to the usual supergravity fields. In the bulk it is related to the metric on AdS^, the metric 
on and the four forms on AdS^ and S^. Following the notations of ||3^, the relevant 
fields sourced by at linear order are 

>AdS rS AdS S ^OOIN 

The fluctuation sources these supergravity fields and these in turn couple to the probe 
D3-brane. To compute the one point function we need to expand the probe D3-brane 
action to linear order in the fluctuations (|3.31| ). The linearized fluctuation contribution 
from the D3-brane DBI action is 

-^jdz dxo dx, d^ [z^h^// + z'h^f + z'htf + h^f + hi,) , 

(3.32) 

where g is the induced metric. 

The linearized fluctuation contribution from the WZ term of the D3-brane is: 

^wz = j a^tpa (3.33) 

As can be seen from the expressions above, our probe D3-brane, which is mainly along 
AdS^, does not couple to linear order to af^^^, the 4-form on . 

The last ingredient we need is the relation between the diagonal fluctuation and 
the fluctuations in the basis ( |3.31|) . This is obtained by solving the linearized Type IIB 
supergravity equations of motion in the presence of the source s = s^Y^ , where is 
an 5*0 (6) spherical harmonic. The solution is given by |]35[] |]36| 

AdS 6 . 4 



h% = 2Asg^p, (3-34) 



AdS ^ _4 yr? 
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where is the usual covariant derivative and V(^Vj,)S denotes the symmetric traceless 
piece of V^VjyS. 

FinaUy we have aU the pieces in place to calculate the correlator. We plug into ( |3.32| ) 
and ( |3.33| ) the expressions for the fluctuations ( p.31| ) in terms of s using ( |3.34| ). For s we 
take the bulk to boundary propagator ( p.30|) with a boundary source 



A fc=-A,-A+2, 



so = Y: E Y^'''(9 = lcp)sP^ (3.35) 



The correlator of a surface operator 0-£ with the Af = 4 SYM CPO OA,k in ( |2.28| ) is given 



by the the derivative of the D3-brane action with respect to the source s^'^{x). It is given 
by 



(Ce) A^cosh^wo 



dz dxo dxi difj 
,A-i^ik((Po-^) (3.36) 



(xq + xf + -2^ — 2dztanh.uo cos{^{j — i/j') + d^)'^+^, 



where we have used that on the D3-brane solution ( |3.9D tp + (j) = (po, where (po is one of the 
parameters characterizing the D3-brane embedding. Integrating over xq and xi and using 
that Td3 = N/2n^ (in the units = 1) yields 

(Ce • OA,k) _ 2^/^VAd^CA,k r ^ ^1 ^A-lgz/c(0o-^) 



dz dtp ■ 



(Ce) TTcosh^Mo Jo Jo (z^ -2dztiinhuocositP-tP') + d^)^+^ ' 

(3.37) 

We make a change of variables from ^ to C deflned by 



2dz tanhuQ 1 
z^+d^ =V ^ 

As z varies between zero and inflnity, C goes from one to cosh uq (at z = d) and back to 
one. The resulting integral is then simplifled by the fact that only terms with a branch-cut 
in C will contribute. This gives 



(Os) 2A/2+%rfA ginh^ Vcosh^ uo - e 



X 



"271 _ ^ikip (3.39) 
X / dtp- 







(C- v^C^cos^)' 



27 



We use now the integral representation of the associated Legendre polynomiallli 



.(A-|fc|)!^|.L 







d^- = 27r^— L-PX^'(C). (3.40) 

(C- A^cosV;) 

After a further change of variables 

a/C^ — 1 = sinh uq sin x , (3-41) 
the correlator ( |3.39|) reduces to the integral 

(A - A;)!v^CA,fce^^('^o-^') ^ 



A!2^/2rf^ 







dx sin"^-^ X V 1 + sinh^ "o sin^ x • (3-42) 



We can now extract the correlators of the surface operator Oj: with any CPO OA,k- 
We write down explicitly the result for A = 2 and 3 which we also did in the gauge 
theory in Section 2 and do in the "bubbling" geometry description in the next section. 
Using the spherical harmonics in Appendix A we have that 6*2,0 = l/v^? C'2,±2 = 1/2, 
Cs,±i = "\/3/4 and 6*3^^3 = 1/2^/2 while the relevant associated Legendre polynomials 
are given by P^{x) = (3^^ - l)/2, P^{x) = 3{x^ - 1), ^3^(0;) = 3Vx^ - l{5x^ - l)/2 and 
P3 (x) = 15(x^ — l)'^/2. Performing the integral in ( p.42|) we conclude that the correlator 
of 0-£ with the CPO's (|2.29D with A = 2 and 3 are given in the probe approximation by 



(02,0 -^e) _ 2 COsh^^o (02,2 -Oe) _ 1 Sinh^^o 2^(0o-^') 

(0e) rf2 ' (0s) ^ d2 e 

(03,1 • Oj:) ^ J_ cosh^ Up sinh^o ,(^,,_^') (03,3 • Oj:) ^ J_ sinh^ up 3z(</>o-V') 

(0e) V2 ' (0e) V3 

(3.43) 

The remaining correlators can be obtained by complex conjugating ( |3.43|) . 

Recall that the probe calculation is done with a single D3-brane and is dual to the 
surface operator O-^ in the gauge theory where the gauge group is broken to t/(l) x SU {N~ 
1). In order to compare to the gauge theory result ( p.30| ) we take A''! = 1 and A2 = A — 1, 



we use the identification (|3.10|) to relate 6*"^° sinhuo to /3i + ^71, and take /32 = 72 = 0. 
Then with z = de'^'^ , we find that (|3.43| ) agrees precisely with (|2.30|) for 02,2 and 03, ±3. 
The other correlators in ( |3.43| ) exactly reproduce the semiclassical result in ( p.30| ) but 



11 We use a slightly modified definition, P^^^x) = ^^{x^ - l)l'='/^(d^+l'=l(x2 - 1)^ /d'^+^^^x) 
which is single valued for argument greater than one. 
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deviate from the classical result of the gauge theory by a quantum correction, proportional 
to the 't Hooft coupling A. The general form of the correlator (|3.42D for Oa./c computed 



in the probe approximation indicates that the quantum corrections to the semiclassical 
gauge theory computation truncate to order X(^~\'^\')/'^ ^ suggesting that the gauge theory 
description of these correlators may be captured exactly by a reduced matrix model, which 
has a finite number of loop corrections. 

3.4- Correlator with Wilson Loop Operator 

We wish to examine the correlator of the surface operator Oj: with the Wilson loop 
operator (|2.38| ). We employ the target-space metric ( p.l2| ) where the brane is at fixed 



u = uq. The Wilson loop is described, as usual, by a semiclassical string surface 
ending on the boundary of AdS^ x along a path related to the specific choice of loop. 
The specific loops we are considering are those studied in [|15|[^, where the relevant 



boundary conditions are further explained. The string worldsheet is not at fixed u, and it 
varies from oo to w = wq, where it ends on the D3-brane probe. In the u ^ (yo region the 
string worldsheet ends in the AdS^ x 5"^ boundary along the desired loop, which is a circle 
along the ip direction at fixedS p. At the same time the string also wraps the direction (p 
along a circle of 5"^ at ^ = 9o in the metric (3T). The extra parameter i/'o which appeared 



in the definition of the loop is the relative phase between the two angles, so that at the 
boundary if; + cp = i/jq, matching the phase of the Wilson loop operator in ( p.38|) . 

Finding the string solution for this configuration is by now a standard procedure. 
With the worldsheet coordinates a and r one assumes the rotationally-symmetric ansatz 

u = u{a), 'il; = T + a{a), 9 = 9{a) , (j, = -t - ({a) , (3.44) 

and one can consistently keep all the other fields constant. The boundary conditions at 
M — s> oo on '0 and (p are such that a{0) — ({0) = i/jq, so that at the boundary we get 
t/j + (p = ipo, which is the phase corresponding to the Wilson loop (|2.38|) . These extra two 
degrees of freedom are required in order to allow the relative phase of ip and (p to change 
as the string moves in the u direction. At the other end of the string worldsheet, where 
the string ends on the probe D3-brane, we have that ip + (p = (po (|3.8|), so a — ( = (po there. 



Due to conformal invariance the radius of the loop in the gauge theory, and the value of p 
here are inconsequential 
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The string Lagrangian in the conformal gauge is 

Two conserved quantities are the momenta conjugate to a and ( 

Pq. = sinh^ M a' , P(^ = sm'^9('. (3.46) 

There are two Virasoro constraints 

sinh^ ua + sin^ ^ C' = . u'^ - sinh^ u{l-a'^) + 9'^ - sin^ ^ (1 - C^) = . (3.47) 

The first one sets = —Pc- By solving the equations of motion for u and 9 one finds 
expressions similar to the second Virasoro constraint for each independently 

2 2 

u'^ - sinh^ u + = a" , _ gi^2 Q ^ ^ _^2 _ ^3_4g) 

sinh u sin 9 

These equations can be solved in general in terms of elliptic integrals, which we will not 
do here, since we will use in Section 4.4 a different coordinate system where the solution 
can be found rather easily. 

Instead, let us focus for now on the simple case with Pq, = —p(^ = which is the case 
when the Wilson loop (|2.38| ) and surface operator 0-£ are in phase, i.e. Vo = <?^'o- One 
solution is for = 1 which corresponds to the Wilson loop on a great circle on 5"^. 

1 TT 

ijj — tpo = —(j) = T , coshw = , 9 = — . (3.49) 

cos a 2 

The coordinate u varies between infinity, the boundary of spacetime, and uq where it ends 
on the D3-brane. After subtracting the divergence from infinity, the finite part of the 
string action is 

S = \ da{u''^ + cosh^ u) = — V^sinhwo . (3.50) 

On the D3-brane worldvolume we have a non-trivial holonomy for the worldvolume 
gauge field = a ( |3.10| ). Since the string ends on the D3-brane along a curve wrapping 
the tj) direction there is a contribution from a boundary term on the worldsheet. It is given 
by 

^ij dip ^ g27ria gj^-j 
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The correlator between a surface operator Oj^ and the Wilson loop is given in the 
semiclassical probe approximation by string worldsheet action. Combining the bulk and 
boundary term and using the relation ( |3.1C1| ) between sinh-uo s-nd the parameters of the 
surface operator Oj: in the gauge theory, the expectation value of the Wilson loop in the 
presence of the surface operator O-e is given by 

■ Oy)_ ^ ^27r|/3 + + 2ma) . (3.52) 



{Oe) 

This probe D3-brane computation reproduces the result of the gauge theory computation 
( Fill ) for gauge group U{1) x SU{N - 1). 



Note that there is another string solution with these boundary conditions |2^, which 
has a = and 

sinh-u = — - — , sin6' = — - — . (3.53) 
smh a cosh a 

This solution does not end on the D3-brane, but is still a solution with the prescribed 

boundary conditions. This solution has zero action. 

Using the coordinates in the next section and the explicit solution in Section 4.4, we 

can evaluate the action for arbitrary tpQ — (^o and 9q. There are two disconnected solutions. 



found already in [|T5|] with action 

cos6'o. (3.54) 

Clearly the solution with a negative sign dominates, but the other one, which is unstable, 
can also be found in the asymptotic expansion of the Gaussian matrix model and we 
therefore retain it. 

The connected solution, with the contribution of the holonomy on the probe D3-brane 
has action 

S = —\/X sin 6*0 cos('i/'o ~ 't'o) sinhwo — 27ria . (3.55) 



Summing over the three saddle points and using ( p.lOj) we get the result 

exp (27rsin6'o(/3cos'i/'o + Tsin-i/'o) + 27ria) + 



(Oe) 

+ exp cos Oq^ + exp cos 6*0 



(3.56) 



This agrees precisely with the gauge theory prediction (|2.41|) for A'^i = 1 , N2 = N — 1 and 
= 72 = in the large limit, including the last term, corresponding to an unstable 
solution. 
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In string theory the difference between a Wilson loop and an 't Hooft loop is simply 
the replacement of a fundamental string by a D-string. Geometrically the surfaces will 
look the same with the only difference being a relative factor oi I/qs = An/gyM the 
on-shell action (|3.50| ) . In this case the D-strings couples through its boundary to the dual 
gauge field on the probe D3-brane, which is now proportional to ry ( |3.1C1| ). Therefore the 
final answer of the calculation of the correlator between a surface operator O-^ and an 't 
Hooft operator in the probe approximation is 

^ =exp ( "2 — sin6'o(/3cos'i/'o + Tsin-i/'o) + luia ) + 

[AtxN ^ \ f AtiN ^\ 

+ exp I cos c/q I + exp I — cos 6*0 I . 

This calculation done with a Dl-brane is exactly the S'-dual of (|3.56| ) accounting also for 
the S duality transformation ct — > ry. 



4. "Bubbling" Surface Operators in A/" = 4 SYM 

In this section we perform calculations with surface operators in A/" = 4 SYM using 
the Type IIB supergravity solutions proposed in as the gravitational description of 
surface operators. We begin with a summary of the main features of the supergravity 
solutions. We then review how to extract correlation functions of gauge theory operators 
from an asymptotically AdS^ x bulk solution and summarize the main formulas. We 
then calculate using the supergravity solutions the correlators of a surface operator Oy, 
with local operators, with a Wilson and an 't Hooft loop. 

4-1. Surface Operator as "Bubbling" Supergravity Solution 

In [|^ the gravitational description of all the maximally supersymmetric operators 
Oy, in a/" = 4 SYM was found in terms of smooth ten dimensional solutions of Type IIB 
supergravity which are asymptotically AdS^ x . These supergravity solutions capture 
the complete backreaction of the configuration of D3 branes in AdS^ x considered in 
the previous section, which describe a surface operator in the probe approximation. 
Here we summarize the main ingredients of the solutions needed for the holographic com- 
putation of the correlation functions. See for more details about the solution and the 
identification with the maximally supersymmetric surface operators in A/" = 4 SYM. 
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By analyzing the symmetries of the maximaUy supersymmetric surface operators in 
A/" = 4 SYM one may observe ^ that the corresponding ansatz for the dual supergravity 
solution can be found by performing a double analytic continuation of the LLM ansatz 
[P] dual to the maximally supersymmetric local operators in A/" = 4 SYM. The Einstein 
frame metric describing a surface operator in A/" = 4 SYM is given by 



ds = y 



2z-l 



2?/ 



V4z^ - 1 



y/4z'^ — 1 

{dx + V)'^ -\ {dy'^ + dxidxi 



22/ 



while the RR five-form field strength is given by 



(4.1) 



y'^T^idx + v) 



2z-l 

y 



-y^ *3d 






y2 



\dx + V) 



2z + l 



y^ *3 d 



z — 



1 -\ 



A d VolAdSs 
A dQs. 



(4.2) 



2; is a function on the space X parametrized by xi, x^ and y, with ds\ = dy'^ + dxidxi and 
y > 0. y is a one-form satisfying dV = *x dz, and hence, the metric and the five-form 
are completely determined by z(xi, X2, y). 

This supergravity ansatz describes a surface operator 0-£ in TV^ = 4 on AdSs x S"^, as 
the metric in the conformal boundary of (|4.1| ) is that of AdSs x S^. The choice of AdSs 
metric determines the surface E on which the surface operator Oj: is supported: AdSs 
in Poincare coordinate corresponds to a surface operator Oy. supported on E = i?^ while 
global AdSs corresponds to a surface operator 0-£ supported on E =S^. 

A non-trivial solution to the equations of motion is obtained by specifying a configura- 
tion of M point-like particles in X . The data from which the solution is determined is the 
particle positions {xi,yi) in X (see Figure 1). Given a particle distribution, the function 
z{xi^X2^y) is found by solving the following differential equation^ 



— — ^ — — j = Yy^'^yi^iv 



yi)5^''\x-xi). (4.3) 



The "charge" Qi = 27ry; associated to these particles is fixed such that the x circle shrinks 
smoothly at y = yi. 
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a) h) 



Fig. 1: a) The metric and five- form flux is determined once the position of the 
particles in X - labeled by coordinates {xi,yi) where y > - is given. The Z-th 
particle is associated with a point Pi G X. b) The configuration corresponding to 
the AdS^ X vacuum. 



This equation is equivalent to the six dimensional Laplace equation with 5*0(4) in- 
variant sources for the variable 

$ = 4. (4.4) 
The solution to (^]^) giving rise to a non-singular smooth metric is given by 

1 ^ 

z{xi,X2,y) = - + ^zi{xi,X2,y), (4.5) 

1=1 

where 

ziix^ X2 v) = ix-xi)'+y' + yf 1 

2^{{x-xir + y^ + yfr-4yfy^ 2" 

The one form V = Vjdx^ can be determined up to an exact form by solving dV = - *x dz 



M 



Vj = -eij V jx'^ -xf){{x-xiy + y^ -yf) 

^ " 2{x - V((^ -xi)^ + y'^ + yfY - Ayfy"^ ' 

Therefore, the metric and five-form field strength are determined by the integer M and by 
the particle locations {xi,yi) for / = 1, ■ ■ ■ , M. 

As explained in 0, further data needs to be given in order to fully specify a Type 
IIB supergravity solution. These metrics have M non-trivial two-cycles Di and a solution 
is fully specified only after we fix the periods of the NS-NS and RR two-form gauge fields 
on the two cycles 

^ ^« / = 1,---,M. (4.8) 



D, 27r 27r 
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y 
Pi 



PA 





Fig. 2: A disk D can be constructed by fibering over an interval connecting 
the "charge" at the point Pi ^ X with {xi,yi) coordinates and the boundary of 
AdSs X 5^ 



Therefore, these "bubbhng" supergravity solutions depend on 4M real parameters 
and there is a one to one correspondence between these supergravity solutions and the 
maximally supersymmetric surface operators in A/" = 4 SYM, which also depend on 4M 
parameters («/, 7^, The explicit mapping is given by 



Bns 
D, 27r 



(4.9) 



2^^ 

The ranks of the subgroups of the gauge group L = Yli^i U{Ni) left unbroken by is 
encoded in the positions in the y axis of the M particles via A^; = , as yf determines the 
amount of five- form flux through the corresponding 5"^. In the units where the AdS^ x 
radius of curvature L = 1, we have that yf = Ni/N. 

These solutions are regular but become singular precisely when the path integral 
describing a surface operator Oj: becomes singular 0. This occurs when the unbroken 
gauge group near E is enhanced, which corresponds in the supergravity solution to colliding 
the charges in Figure 1. When two charges collide an S'^ in the geometry shrinks to zero 
size and the solution is singular. 
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4-2. Gauge Theory Correlation Functions from Ten Dimensional ''Bubbling" Solutions 

Here we summarize the main ideas and formulas involved in extracting the correla- 
tion function of local operators in the boundary gauge theory from the ten dimensional 
"bubbling" supergravity solution. We refer to for more details and derivations. 



The asymptotically AdS^ x 5"^ solutions of Type IIB supergravity describing maximally 
supersymmetric surface operators 0-£ in A/" = 4 SYM are ten dimensional. Holography is 
on the other hand an equivalence between quantum gravity with AdS^ boundary condi- 
tions and a four dimensional field theory on the boundary Therefore, given 



a ten dimensional asymptotically AdS^ x solution of Type HB supergravity, one must 
systematically reduce the ten dimensional solution to a solution of five dimensional grav- 
ity coupled to — in general — an infinite number of five dimensional fields. Once this 
reduction and a local five dimensional action has been found, one can apply the usual 
holographic rules |^ |Q and use the bulk five dimensional description to compute the 
correlation function of boundary operators. In order to get the renormalized correlation 
functions in the gauge theory one must appropriately renormalize the bulk gravitational 
action using holographic renormalization (for a review and more references see 



The first step in constructing the five dimensional gravity action from which one 
can holographically compute gauge theory correlation functions is to decompose the ten 
dimensional solution into AdS^ x plus the deviation of the solution away from AdS^ x S^. 
One then expands the deviation in a complete basis of SO{6) spherical harmonics which 
yields an infinite collection of five dimensional fiuctuation fields. The decomposition into 
harmonics of the the metric and five-form deviations is given by 

Ka{x.y) = Y,iB[l^^x)Yt{y) + B[l^^{x)D^Y'-{y)) ^ 
hab){x.y) = Y.{^[l^{x)Yl^^^{y)+ct>'(;:,^{x)D^,Yif{y) + 

K{x,y) = J2^'^{^)y'Hy), 

and 

fabcd^.ix,y) = y^{D^h\lAx)eab,d^D,Y'^{y) + {K'^ - A)h'^- {x)eabcd^Yt {y)) 

^ > ^ (4.11) 

fabcdeix.y) = ^b^^l^{x)A^^eabcdeY^^{y), 
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where are AdS^ coordinates while are coordinates. Y^^, Y^^ and Y^^^^ are 
scalar, vector and symmetric traceless SO{Q) spherical harmonics respectively while — A''^^ 
and — are the eigenvalues of the d'Alembertian on scalar and vector harmonics (see 
Appendix A for more details on the relevant spherical harmonics). One must then solve the 
ten dimensional Type IIB supergravity equations of motion perturbatively in the number 
of fluctuations. 

In [^5[ the linearized analysis was performed in the de Bonder gauge which retains 
the physical fluctuations. For our solutions it is convenient to work with gauge invariant 
variables as our "bubbling" solution is not in the de Bonder gauge. The perturbative 
analysis using gauge invariant variables has been carried out by Skenderis and Taylor in 
| ]4(j| j, where they showed that the equations of motion for the gauge invariant variables 
indeed coincide with the equations of motion in the de Bonder gauge if the fluctuations 
in the de Bonder gauge are replaced by their gauge invariant counterparts. The gauge 
invariant variables to linear order in the fluctuations are given in terms of the fluctuations 
in (TOdiTTDby i9| 













= bi\ - 








= b!' 








= 6^5 - 

At 


1 




2(A^5 - 


4) 




= 6^'\ 







(4.12) 



The gauge invariant fluctuations that diagonalize the equations of motion to leading order 
in fluctuations and that couple to the CPO's 0^ with A = 2,3, the SO{Q) /^-symmetry 
current and the stress-energy tensor Tmn of A/" = 4 SYM are 

80 V ) 



100 V I (4.13) 

"'^1L' "-/it/ ^y^i,'^ ■ 



In order to calculate correlation functions in the boundary gauge theory from the bulk 
description one must flnd a local flve dimensional action that controls the dynamics of the 
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gauge invariant fluctuations dual to the gauge theory operators. This requires flnding a 
non-hnear map between local flve dimension flelds and the gauge invariant fluctuations 
obtained by expanding the ten dimensional equations of motion. This mapping is con- 
structed such that a solution of the flve dimensional equations of motion yields a solution 
of the ten dimensional Type IIB supergravity equations of motion. This mapping can be 
constructed perturbatively in the number of fluctuations. This non-linear map between 



local flelds and gauge invariant fluctuations has been constructed in . 

Once a local flve dimensional action yielding the appropriate flve dimensional equa- 
tions of motion is constructed, the correlation functions of renormalized boundary op- 
erators can be computed using holographic renormalization. In order to construct the 
renormalized correlation function of boundary operators the bulk action must be enriched 
with a set of boundary terms to make the on-shell bulk action flnite. Despite the fact that 
the flve dimensional action in general contains an inflnite number of flelds, only a flnite 
number of terms in the action are infrared divergent. It is only the bulk flelds dual to 
low dimension operators in the gauge theory that require boundary counterterms. Once 
a complete set of counterterms for the flve dimensional action are constructed, one can 
calculate the correlation functions of boundary operators by taking functional derivatives 
of the renormalized bulk gravitational action. 

Here we are interested in computing the correlation function of the CPO's 0^ with 
A = 2, 3, of the /^-symmetry current and of the stress-energy tensor Tmn with a surface 
operator O-^ in A/" = 4 SYM. The bulk action that captures the degrees of freedom dual 
to these operators is given by 11110 



T\j2 r / /? 1 1 



262 

2 



3 16 ^'^ ^'^ 



(4.14) 



where F^^, is the fleld-strength of a vector potential and G^ii/, O, and couple to 
Tmm C^A=2' ^A=3 ^rn respectively and a is an SO{Q) vector index. The boundary 
counterterms required to have a flnite on-shell action and from which the correlation func- 
tion of renormalized operators can be computed can be found by solving the equations of 
motion of ( [4.14| ) near the boundary of AdS^ . They are written down in . 



This action is written, as throughout the rest of the paper, in units where the radius of 
curvature of AdS^ x is L'^ = AirgsNlt = 1. 
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In order to calculate the one point function of these local operators using the bulk 
description it is convenient to write the five dimensional metric G^i, in the Fefferman- 
Graham form 

dU"^ 1 

dsl = -^ + j^ {G^O)mn + U^G(2)mn + " " ■) rfx'^dx", (4.15) 

where in our case G(-o) is the AdS^ x metric on the conformal boundary at t/ = 0, where 
A/" = 4 SYM lives. The bulk fields O, \E' and also have a near boundary expansion. It 
is given bylll 

^ = s^\x,U) = U^ [s^^]^Y^ + ... (4.16) 

A«(x,t/) = t/2 [a;]^{x)y^- + ..., 

where [ ■ ] is the normalizable mode of the field. 

Evaluation of the on-shell one point function of the renormalized bulk action shows 
that the vacuum expectation value of these operators is captured by the normalizable mode 
of the dual bulk field. In order to calculate the one point function of a unit normalized 
operator, we must divide the one point function obtained from supergravity by the square 
root of the two-point function obtained from supergravity , such that the supergravity 



computation of the vacuum expectation value corresponds to that of a unit normalized 
operator with two-point function given by ( |2.27| ). After taking this into account, we have 
that the one point function of the CPO's (9^ with A = 2, 3 and of the /^-symmetry current 
are given by 

(^30o.=|^[^"]3^ (4-17) 

(r) - - ^ U«l 

- 247r2 ™2 • 

The computation of the vacuum expectation value of the stress-energy tensor Tmn is 
more subtle as the bulk scalar field Q{x, U) enters in the renormalization procedure for the 



In our solutions the bulk fields do not have a source term and therefore we do not write it. 
In a general solution a non-normalizable mode has to be added to these asymptotic solutions. 
The current is normalized so that it has the canonical two-point function |43]. 
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five dimensional metric. The formula for the vacuum expectation value of the stress-energy 
tensor is given by MB 



{Tmn)o^ - ^ (^fi'(4)mn " g [s^^] 3 i^"^^] 2 9{0)mn 

(4.18) 

In this formula gi^.)rnn denotes the various terms in the near boundary expansion of the ten 
dimensional metric in the Graham- Fefferman coordinate z in ( |4.15 ). 



As explained in Section 2, the expectation value of a local operator {0)q^ in a non- 
trivial state of the gauge theory on AdS^ x S"^ corresponds to the two-point correlation 
function of the surface operator Oy, with the local operator (9 in A/" = 4 SYM on R'^ 

■ 0> (4 19) 

Therefore, we can calculate these correlators computed in the previous sections in other 
regimes by using the supergravity description of surface operators. 

4-3. Correlator of Surface Operators with Local Operators from "Bubbling" Solutions 

We can now compute the correlation function of a surface operator Ot, with the CPO's 
Oj^ of A = 2, 3, with the U{1) current and with the stress-energy tensor Tmn of A/" = 4 
SYM using the holographic correspondence. The analogous calculation for local operators 



has been carried out by Skenderis and Taylor in [^. As summarized in the previous 
subsection, the computation of the correlator of 0-£ with these local operators requires 
extracting the appropriate coefficient of the near boundary expansion of the corresponding 



fluctuation mode in (^4.17|) ( ^T8| ) from the "bubbling" supergravity solution (^J]) 



Before looking at the fluctuations, we write down a reference AdS^ x S^ solution, 
which is the metric (|4.1D(^?^ with a single particle source. If the particle position is at 



{x^'^\ y*-'^-' = 1), then in the coordinates 

-x^°^ +z(x2 =re*" 
r = R sin 9 



y = a/1 + cos 6* (4.20) 



a = ip + (j), 
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we have that: 



(0) _1R^ + cos^ e + 1 



(4.21) 



+ sin^ 9 

Inserting ( |4.21D into ( |4.1D results in the AdS's x S"^ metric with AdS^ fohated by AdS'axS^ 
shces0 



1 + R^)ds\^s-, + XTb^ ^ ^^^^^ ^^^^ ^^"^ 



ds' 



The five-form flux ( [4 .21 ) is given by 



2 j„/.2 



,2 /) J i2 



(4.22) 



Fs = i?(i?^ + l)dR A di/j A d VoUdSa + cos^ 6* sin 9d(j) A dO A dn^. 



(4.23) 



The supergravity solution ( ^72|) describing a maximally supersymmetric surface 
operator Oj: is characterized by a choice of particle distribution (x;, yi) with / = 1, ■ ■ ■ , M. 
This determines the function z{xi,X2i y) and the one form V. For computational purposes 
we flnd it useful to use the function $ ( ^I^) , which is given in terms of z{xi,X2,y) by 



(4.24) 



We then rewrite ^{xi,X2,y) and Vi{xi,X2,y) as deviations from the reference AdS^ x 5"^ 
solution 



$ = $(0) + A$ , ^^^^ = 
V = + AV. 



z^ 

y2 



(4.25) 



It turns out that the most useful reference metric, around which the fluctuations are 
minimal, is the one centered around the center of mass deflned bylll 



M 



M 



(4.26) 



1=1 



1=1 



Furthermore, to simplify the calculations, henceforth we take that x^^^ — 0. 



This is the same as ( p.l2| ) with R = sinhu. 

Recall (|4.9| ) that yf = Ni/N, and is therefore a natural weight for the point xi. 
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Now, since the function A$ solves the six dimensional Laplace equation with SO (4) 
symmetry, all the fluctuations can be expanded in a basis of SO (4) invariant harmonics 
F^''^ of SO {6) (see Appendix A for a summary of their properties). Therefore 

A«=^ t Ai.. j"-^lf:^\ (4.27) 

A>2 fc=-A,-A+2,--- 

where 9,(1) are coordinates on the S^ ( |4.22D . The A = 1 term in ( [4.27| ) vanishes in the 
center of mass frame which we have adopted. 

By doing the calculation in the center of mass frame (gj^) with x^^^ = 0, we get from 
the M particle solution (|]5|)(|]6D thaM 

•2 2/^-l^„.2 



A$2,o = 4v^5^ f 
1=1 ^ 



TT- Vl 



M 



A$2,±2 = 6eT2^V'^(^^^^±,^^2: 
1=1 

M 

A$3,±i = 8V3e^'^ Y.^xi,i ± ixi,2){xf - yf + -)yf 



(4.28) 



1=1 

M 



A$3,±3 = 8v^eT^3^ J2^xi,i ± ixi,2: 



1=1 



Now to calculate the desired expectation values of Oa./c, Jm Tmn we use the rela- 
tion ( [4.171) (|4.18|) between the expectation values and the fluctuations [s^^] ^ , [s^^] 3 , [^Ji] 



2 



and g(.)mn- These in turn are related to the moments of A$ which we have just calculated 
by following the procedure in Section 4.2 as we outline now. 

We need to find the asymptotic form of the metric and five-form and then expand 
it in spherical harmonics as in ( [4.10| )( |17TlD . This is done by plugging the coordinate 
transformation (|4.20| ) into the M particle solution ( [4.5|)(|4.6| ) (|4.7| ) and expanding for large 
R. In the multi-particle case this will not yield a manifestly asymptotically AdS^ x 5"^ 
metric. This can be fixed by modifying the one- form V ( [4.7| ) . Recall that V is determined 
from z by solving the equation ydV — *xdz and hence is only determined up to an exact 



These expression will also be valid away from the center of mass frame, with a subtlety ap- 
pearing only in A$3,±i, where the factor of 1/2 needed to be determined by a separate calculation. 
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form. In order to make manifest the AdS^ x asymptotics we add to V in ( [4.7|) an exact 
one form 

V ^V + duj. (4.29) 

The exphcit expression for uj that yields the desired asymptotics is given in Appendix D. 

Now we can compute the various fluctuation modes. The warp factors in the metric 
( [4.1|) can be parametrized as foUows 



(4.30) 



V4z2 _ 1 2|/ (i22 + sin^5/) 



The functions A^B^C^D encode the deviations of the metric away from AdS^x'^^ . We 
express them in terms of A$, which scales (in the center of mass frame) like 1/ for large 
R i W^ - Up to 0(l/i?^) they are given by 

A = -R^A^ + R^ sin^ 6'A$ - -R^{A<^)^ 
2 8 

B = --R'^A^ - R^ sin^ OA^ + -R^{A^f 

C = --R'^A^ - R^A^ + -i?S(A$)2 
2 8 

D = -i?^A$ + R^A^ - -i?^(A$)2. 
2 8 



We also expand the flve-form flux ( [4.2|) in a large R expansion. The expression for A$ is 
obtained from the explicit expression for the M particle solution ( [4.5|) (|4.6|) . 

We now need to use (^TT^) (^]TT|) to calculate the fluctuation modes that appear in 



( [4.17| )( |lT8D which are needed for the computation of our correlation functions. Recalling 



the expression for the gauge invariant fluctuations ( |4.13| ), we flnally have that 







k = 2,0,-2 






/c = 3, 1,-1,-3 




-^(A*).„. 





(4.32) 



Now we recall the expressions for the moments of A$, which we calculated in ( |4.28| ) 
and can immediately write down the supergravity answer for the correlator of a surface 
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operator Oj: with the CPO's OA,k ( |2.29| ) with A = 2, 3. Writing the answer in gauge 
theory variables using (^4.9| ) we get that 



(0E-02,o) _ 1 8n' (^^^ /,^2 , 2^ A N-Ni 



IVl / 



(0e) v^A ' V 47r2 2iV 

(4.33) 

(^E ■ 03,i) 1 Stt^ ^ / ^ A iV-2iV, 



{f3i + hi] 



{Oj: ■ 03,3) 



1 o 3 



The remaining correlators can be obtained by complex conjugating ( [4 .331 ). 

To compare with the probe brane calculation (|3.43| ) in the previous section, we consider 
the result of the supergravity calculation ( [4.33| ) for A^i = 1 and N2 = N — 1. Then we 
impose the center of mass condition, which is valid for the probe calculation and finally 
take — > 00 limit. The results are in precise agreement. 

We can compare with the gauge theory for arbitrary surface operator parameters 
{Ni, (3i, 'ji). As in the case of the probe calculation ( |3.43| ), the supergravity computation 
matches exactly the gauge theory result (|2.30| ) for operators with A = |A;| and the semi- 
classical gauge theory result in ( |2.30| ) exactly matches the leading term in the supergravity 
computation (|4.33| ) for the rest of the chiral primary operators. The supergravity com- 
putation suggest that the loop corrections to the gauge theory result for the correlator 
between (9^ and OA,k truncate at order Xi^-l'^D/'^. 

The vacuum expectation value of the U{1) current in given byi^l 

The first term agrees precisely with the gauge theory result ( |2.36D and the second captures 
the one loop correction. 



In computing the correlator with we have used that the supergravity solution is U (1) 
invariant and that this U{1) is the diagonal sum of the [/(I) C 50(1,5) generated by which 
rotates the space transverse to S and the U{1) C 50(6) subgroup of the i?-symmetry generated 
by Jt. The [/(I) symmetry implies that (j^)^ = {jt) ^ ■ 
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We now proceed to calculate the correlator of Oj: with the stress-energy tensor Tmn- 
We first calculate the near boundary expansion of the metric along the AdS^ directions. 
It is given by 

V 4v^i?2 6 X 24i?2 ; R^ + l\ 3 X 24/^4 J 

where we sum over the U{1) charge /c = 2, 0, —2. 

The correlator of 0-£ with can be obtained from formula ( |4.18| ). This requires 
bringing the ten dimensional metric into the Fefferman-Graham form ( [4.15| ). In order to 
bring the metric (|4.35| ) into this form we must perform the following change of coordinates 



1 /I (A$2,fcA$2,-fcj 



R V 4i?2 24 X 24 

U \ 4 24 X 24 

In the Fefferman-Graham gauge the metric ( ^4.35| ) is given by 

rf,^ = rf^^ I J_ _ 1 + 18 + 24v/3A$2,o + (A$2,fcA$2,-fe) ^,2^ ^ dU'^ 



(4.36) 



t/2 2 288 / t/2 

,1 1 18-8v^A$2,0 + (A$2,fcA$2,-fc).,2 

+ ^^^^^53 1^+2+ 288 ^ 



(4.37) 



Now we can read off from ( [4.37| ) the various terms in the metric ((7(4), 5^(2) , 5^(0)) that 
enter in the formula of the correlator of Oj: with the stress-energy tensor T^nn ( [4.18| ) . The 
correlator is given by 



^ Ar2 / 3 I \ 

A^2 / 1 1 \ 



(4.38) 



where gab is an AdSs metric and the explicit expression for A$2,o is given in ( [4.28| ). More 
explicitly 



3A^2 2 ^ , o. X N -N, 



2 / 
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(4.39) 



The A independent term encodes the Casimir energy in the gauge theory. The leading 
term in A precisely agrees with the gauge theory calculation in the semiclassical regime 
( p.35|) . The last term can be interpreted as the one loop correction to the gauge theory 
result (|2.35| ), and suggests that the correlator of O-^ with the stress-energy tensor Tmn 
receives only a one loop correction. 

We note that the trace is {T^)q^ = 0, which reproduces the fact that the conformal 
anomaly of TV = 4 SYM on Ad S3 x vanishes. 

4-4 ■ Correlator with Wilson and 't Hooft Loop 

We turn now to calculating the Wilson and 't Hooft loop expectation value in the 
supergravity background describing a surface operator Oj:. As usual in the holographic 
dual, the Wilson loop is described by a classical string [|3^[Q. As in the calculation in 
the probe approximation in Section 3.4, we have to identify the boundary conditions on 
the string stemming from the parameters of the Wilson loop VFeo,^(j (|2.38|) and 't Hooft 
loop Unlike the probe approximation, here the string has no D3-branes to end 

on, rather it will wrap a non-trivial cycle in spacetime and close smoothly on itself, as we 
explain shortly. 

To understand the boundary conditions we look at the asymptotic form of the metric 
( [4.20|)(|4.22|) and impose the same conditions as in the probe calculation in Section 3.4. 
At the boundary, where 00, the string should wrap a circle of at 9 = 9q, so the 

coordinates r and y in ( ^4.20|) should diverge with a fixed ratio r/y tan^o- The string 
wraps circles both in AdS^ and on S^, parametrized by i/j and (p with a relative phase 
a = + 4) = if^o, corresponding to that of the Wilson loop VF^cV'o ( P-38|) . 

To get a smooth string solution with these boundary conditions, it has to have the 
topology of the disc and close smoothly on itself. Recall that the x circle shrinks to zero 
size at each of the sources {xi,yi), so the string should extend from the boundary to one 
of the sources. 

We take the ansatz where only x depends on r 

Xi^Xi{a), y = y(a), x = r + r]i(^). (4.40) 
Using the metric ([4.1j ), the string action in the conformal gauge is 

^ = ^ / {vw=i ('"' + + + ^^^(»" +^-^-^') • f^-"' 
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where V^^ is the puUback of the one-form V in the a direction. 

There are two Virasoro constraints for the diagonal and off-diagonal pieces in the 
stress-energy tensor 

v' + v^ = o 

, , 4,^ (4.42) 

The first equation will give ry, once we solve for V (which depends only on (x, y)). Then 
plugging the second Virasoro constraint into the equations of motion for the xi and y one 
finds 



These equations mean that the ratios of y' , x'l and X2 are constant along the worldsheet. 
Therefore the string is given by a straight line in this space! 

These equations are actually unnecessarily complicated, since they are in the conformal 
gauge. We can switch to the Nambu-Goto formulation, the the string action is very simple 



27V 



j da dr y'"^ + x'^x'^, (4.44) 



In this formulation rj does not appear and can be set to zero. The action is clearly just 
the length functional in (x, y) space, which is solved by straight lines. These worldsheets 
correspond to the disks spanned by the straight lines in Figure 2. 

As mentioned above, the solution will have to extend from the boundary of spacetime 
to one of the sources, where the x circle shrinks away. The length of the ray emanating 
from the source (x^, yi) in the direction of the boundary point 

(i?sin 6*0 cos t/'o, -Rsin6'o sinT/^o, -Rcos6'o), (4.45) 

with R 00 is given by 

R — yi cos 6*0 — sin6'o(xi,i cos -00 + xi^2 sini/jo) . (4.46) 

Recall that the "bubbling" supergravity background also has an NS-NS two-form 
gauge field so in addition to the Nambu-Goto term, the action includes the fiux through 
the worldsheet, which according to ( [4.9| ), for a string ending on the source / is —2nai. 
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After integrating over r and removing the usual infrared divergence, which here is 
proportional to the cutoff R, the total action is 



S = - {yi cos^o + sm9o{xi^i cos -00 + xi^2 sin-i/^o) + ^niai) . 



(4.47) 



Finally we use (^]^) to rewrite the result in terms of the parameters from the gauge theory 
description, and sum over the different classical solutions, one for each source to get 



M 

^exp 

1=1 



N 



cos 6*0 + 27r sin 6*0 (A cos '00+7 sin ifjo) + 2mai 



(4.48) 



As was mentioned in Section 2.4, when expanding the exact result of the matrix model 
that captures the supersymmetric Wilson loops one finds two saddle points, one with a 
positive exponent of \/~X and one with a negative one. In this was reproduced from a 
second string solution, an unstable one. Such solutions exist also in our case. The reason 
is that the (x, y) space is restricted to y > 0, but the plane at y = is not a boundary. In 
fact, if we consider a string with the ansatz ( [4.40|) , it would simply be reflected from the 
y = plane. 

We therefore have to amend ( [4.46| ) . There isn't a unique geodesic from the boundary 
point ( |4.45| ) to each of the sources {xi, yi), but two. The second one is reflected through 
the y = plane and its divergent length 



R + yi cos 6*0 — sin6'o(x;^i cosifjo + xi^2 sini/'o) • 



(4.49) 



In addition to ( [4.48| ) we therefore find another sum coming from the other saddle points, 
which should all be unstable 



M 

^exp 



1=1 



N 



cos 6*0 + 2tv sin 6*0 (/9z cos 'i/'o + 7 sini/'o) + 2Tviai 



(4.50) 

Note that the sign change is only on the first term in the exponent. 

The sum of ( ^4.48| ) and (|4.50| ) is in precise agreement with the result of the gauge theory 
calculation (|2.41D including our expectation for the term coming from self-contractions of 
the extra scalar cp^. We could not reproduce the prefactors which would require a string 
one-loop calculation to determine. 
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This is immediate to derive by using the method of images with a source at the point (xi , —yi). 
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For an 't Hooft loop Tg^^^^ we get the same action, except that the string tension is 
multiphed by a factor oil/ Qs = An / Qym instead of couphng to the NS-NS field it gets 
a term from the RR field. The correlator is 



M 

^exp 

1=1 



4:71x1 — COS 6*0 H — 9 — {Pi COS '00 + 7 sin ipo ) sin 6*0 + 27Tirji 



A 



9ym 
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to which again we should add the contributions of the unstable saddle points. This result 
is 5'-dual to ( |4.48|) , accounting for ai rji and also matches the general expectations for 
the behavior of the 't Hooft loop as discussed in Section 2. 



5. Discussion 

In this paper we have began to study the properties of surface operators 0-£ in A/" = 4 
SYM performing various explicit calculations with them. We have studied the properties 
of surface operators 0-£ using three different realizations, consisting of the gauge theory de- 
scription, the probe D3-brane description and the "bubbling" supergravity description. We 
have found that a rather nice story emerges between these different realizations, allowing 
for precise matching between quantities computed using different descriptions. 

For the expectation value of the surface operator with E = i?^ we found from the 
gauge theory and from the probe approximation that the expectation value is unity. It 
would be interesting to learn how to perform this calculation in the framework of the 
"bubbling" solutions, but it would probably be more advantageous to do it in the case of 
the surface operators in the six dimensional theory dual to M-theory on AdSj x S"^ |57], 
since unlike our case, those operators do suffer from conformal anomalies. 

For the correlator of a surface operator Oj: with a chiral primary operator (9A,fc and 
with a Wilson loop VFgo.V'o calculation in the three different realizations of 

the surface operator Oj: and found remarkable agreement. The leading behavior agreed 
between all three regimes. Particularly exciting is the case of the Wilson loop, whose expec- 
tation value depends in general on the parameters of the surface operator (A^^, 7^, ai), 
those in the Wilson loop (6*0, i^o) and the 't Hooft coupling A. 

In all these cases we restricted to a purely classical calculation on the gauge theory side, 
but we expect that those correlators will receive some quantum corrections. Remarkably, 
the computations performed in the probe approximation ( |3.43| ) and using the "bubbling" 
supergravity description precisely reproduce the semiclassical gauge theory result ( |2.30|) . 
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Moreover, we have found that the correlator of Oj: with some local operators have ex- 
tra terms in ( |3.43|) ( [4 .331) which are of higher order in the 't Hooft coupling compared to 
the semiclassical gauge theory result ( p.30| ), and that nicely organize into a perturbative 
expansion in A. Given these results, our expectation is that there are perturbative cor- 
rections to ( |2.3U| ) which, nevertheless, terminate at at finite order in the expansion in A. 
In particular, our results suggest that the correlator of Oy: with a chiral primary operator 
Oa k terminate at order reproducing (|3.43| ) ( ^733|) . It would be interesting to 

reproduce this structure of the corrections directly from the path integral description of 
surface operators 0-£. 

In the case of the Wilson loop VFeo,V'o appearance of terms of the form 

±^XNi/N ( CT) , which we attribute to the extra scalar (p^ in (|2.38| ). Again, we did 
not pursue the exact perturbative calculation of its effect, but we anticipate that it will 
give some modification of the matrix model of and [^, perhaps along the lines of 



Hiiiii. 

There exists a definition for the gauge theory in the presence of the surface operator 
Oj: also beyond the classical level. The path integral has to be done in the presence 
of the prescribed singularity. But concrete calculations have never been performed with 
that prescription (prior to this paper no calculations have been done even in the classical 
limit). Our results raise the prospects of being able to do such calculations in practice and 
in the supersymmetric cases considered in this paper possibly improving on our results and 
computing all quantum corrections. It would be interesting to use localization techniques 
to study whether the path integral of A/" = 4 SYM in the presence of a surface operator 



becomes that of a reduced model, extending the recent result of Pestun , who has 
derived the matrix model in [^||2^ by analyzing the Af = 4 SYM path integral in the 
presence of a supersymmetric circular Wilson loop. 

The fact that a supersymmetric surface operator O^, in A/" = 4 SYM has multiple 
realizations where calculations can be done with precise agreement has an analog in the 
case of supersymmetric Wilson loops. These can be described in a variety of ways, pertur- 
batively in TV = 4 SYM [|4| 



as strings in AdS^ ||4^ [jT^ , as a configuration of 

D3-branes ||5^[^1|[0 or as a configuration of D5-branes [^||5^], and finally as asymptot- 
ically AdS^ X "bubbling" supergravity backgrounds . This paper aspires to 



Note that A — |fc| is always even, so we get a standard perturbation expansion in the 't Hooft 
coupling A. 
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be a first step in generalizing the impressive series of comparisons done for Wilson loops 
between gauge theory calculations and probe calculations. Much dynamical information 
can also be extracted from the "bubbling" geometries for Wilson loops [|5^ , extending 
the results obtained for Wilson loops in the context of topological string theory [E9[ using 



"bubbhng" Calabi-Yau's [30||61 



It would be interesting to perform similar computations to the ones in this paper for 
the surface operators of order type constructed in 0]. One can also consider maximally 
supersymmetric domain walls in A/" = 4 SYM. For these the gauge theory realization is 
given in terms of the low energy limit of the Hanany-Witten like setup involving D3, D5 and 
NS5 branes |^2| , the probe brane description is given by the branes in AdSr, x found by 
Karch and Randall and the "bubbling" description is given by the smooth Osp(4|4) 
invariant supergravity solutions found in ||63[^[^. The gauge theory description of 
supersymmetric domain walls in A/" = 4 has been considered in , and a recent 

systematic study of supersymmetric boundary conditions has appeared in |^ . 



The maximally supersymmetric operators in A/" = 4 SYM are of great interest. Unlike 
the vacuum, there are many of them, with their classification still not complete. Yet, they 
posses a large degree of super symmetry, sixteen supercharges, and therefore their dynamics 
are under good control, in some cases to all orders in perturbation theory. 

The maximally supersymmetric local operators in A/" = 4 SYM are captured by a 
reduction from the four dimensional gauge theory to a matrix quantum mechanics of a 



single Hermitean matrix |j7T|] , and are classified in terms of Young tableau. This matrix 
model can then be solved in terms of free fermions whose phase space has been identified 
by LLM in terms of boundary conditions for the "bubbling" supergravity description . 

The maximally supersymmetric Wilson loops in A/" = 4 SYM are also of great in- 
terest and unlike the case of local supersymmetric operators, circular Wilson loops get 
contributions from all orders in perturbation theory. In fact the circular Wilson loop in 
A/" = 4 SYM is captured by a zero-dimensional matrix model |2^]|j25|, a result that has 
been derived recently using localization |^^. This matrix model reproduces correctly the 
expectation value of the loop as calculated by strings, D-branes ||5^ |]72| and also repro- 
duces properties |j73] of the "bubbling" supergravity description |]5^ . Modifications of 



the matrix model also encode the interaction of Wilson loops with local operators [p8 |. 

It is tempting therefore to try to see if there exists a matrix model that captures the 
essential properties of the maximally supersymmetric surface operators. We speculate in 
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the following paragraphs over a possible candidate. It would be interesting to use localiza- 
tion to derive the reduced model capturing the properties of maximally supersymmetric 
surface operators in A/" = 4 SYM. 

Let us stress from the onset that it is not completely clear what properties this ma- 
trix model should capture. In the case of the supersymmetric circular Wilson loop, it 
calculates its expectation value, which may translate to some conformal anomaly in our 
case. Alas we saw that the maximally supersymmetric surface operators 0-£ do not have a 
conformal anomaly. But one property that we do expect it to capture is the classification 
of surface operators in terms of a Levi group L = U{Ni) C U{N) and the parameters 

{ai, (3i, 7z, r}i) with / = 1, ■ ■ ■ , M. 

One can consider a zero- dimensional complex matrix model with only a commutator- 
squared action 

I?$exp ( ^— Tr[$, $]2 J . (5.1) 
\9ym J 

A vacuum of this matrix model is given by a matrix that commute with its complex 
conjugate. Such a matrix can then be diagonalized with eigenvalues /3; -|- 27^, and the 
remaining symmetry, which depends on the degeneracy of those eigenvalues, would be 
exactly such a Levi group L = U{Ni). 

This matrix model has an interesting feature. Due to the lack of a quadratic term 
in the action, when the path integral (|5.1| ) is written in terms of its eigenvalues including 
the proper gauge fixing, the usual Van der Monde determinant is exactly cancelled by a 
similar contribution from the Fadeev- Popov determinant []7^ [|73[ . 

Normally the eigenvalues are repelled from each-other due to the Van der Monde 
determinant, and this repulsion can then be captured by assigning the eigenvalues Fermi 
statistics. For our matrix model, there is no such eigenvalue repulsion and the vacua may 
well consist of a large number of coincident eigenvalues. 

It is therefore natural to assign to the eigenvalues of the matrix model (|5.1| ) Bose 
statistics, which nicely captures the structure of boundary conditions of the "bubbling" 
solutions for surface operators Oj^. The structure of the "bubbling" supergravity solu- 
tions for the surface operators 0-£ is similar to that of the maximally supersymmetric 
local operators in A/" = 4, where there is also a 3-dimensional space X, with coordinates 
{xi, X2, y)- The main difference between the "bubbling" solutions for local operators and 
surface operators is, however, in the boundary conditions for the Laplace equation ( [4.3| ), 
which determines the metric and five-form. 
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In the case of local operators the supergravity data which determines the solution 
is captured by binary information on the y = plane, which is interpreted as the 
occupation number of the free fermions. In our case, corresponding to maximally super- 
symmetric surface operators 0-£, the data which determines the solution is given by the 
position of point-like sources away from y = (see Figure 1). Projecting this picture on 
the y = plane, it can be interpreted as a distribution of bosons. The distance in the y 
direction which represents the number of coincident D3-branes (or eigenvalues) through 
yf = Ni/N, is interpreted as their occupation number and is indeed quantized in terms 
of positive integers. The "bubbling" solution for local operators is encoded in terms of 
fermionic eigenvalues, while the "bubbling" solution for surface operators is encoded in 
terms of bosonic eigenvalues, a property captured by the matrix model (|5.1|). 

One can enrich the matrix model (|5.1| ) by inserting other observables in the integrand 
made out of the fields $ and such as chiral primary operators. At least at the semi- 
classical level this matrix model captures the correlators of the surface operator Oj: with 
chiral primary operators. 

While this picture is satisfyingly consistent, the matrix model ( |5.1| ) has some short- 
comings in describing the surface operators Oy.. For example, it does not capture all the 
details of the vacuum structure. In particular, it does not capture the parameters ai and 
rji needed to fully characterize a surface operator Oe, and that enter in the correlators od 
O-E with the Wilson and 't Hooft operators. 

One may be even more ambitious and hope to find a reduced model that captures the 
interactions between the different objects considered in this paper, including Wilson and 
't Hooft loops. An example to keep in mind is the normal matrix model, which in [p8 



was shown to reproduce the interaction of maximally supersymmetric local operators and 
maximally supersymmetric Wilson loops. 

A possible avenue is to replace the 0-dimensional matrix model ( |5.1| ) with matrix 
quantum mechanics, as is the case for the local operators. Given that instead of the 
theory on S'^ x R we are dealing here with A/" = 4 SYM om AdSs x S'^, we propose to 
Wick-rotate and compactify the time direction. We propose to consider matrix quantum 
mechanics of a single complex matrix at finite temperature 



Tr 



(5.2) 
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Here the direction is compact and the covariant derivative is -D^$ = 9^$ — $]. In 
a vacuum of this theory, one can choose a gauge where is diagonal, and it would then 
be characterized by the parameters ai in (|2.1j ) that label a surface operator Oy,- 

In this matrix model it is easy to define, in addition to the local operators, a Wilson 
loop operator wrapping the tj) direction, which would be similar to (|2.38|) , at least for 
^0 = 7r/2. More degrees of freedom would be required in order to capture the dynamics of 
the extra scalar field cf)^ . We leave these speculations for future investigation. 
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Appendix A. Spherical Harmonics 

In this Appendix we summarize some of the properties of SO (6) spherical harmonics 



(A.l) 



(see for more details). The scalar and vector SO{6) spherical harmonics satisfy the 
following equations, 

Ay^i=A^iy^i, A^i=-A(A + 4), A = 0, 1,2, ■■■ 

= y/s ^ f^h = _(^^2 ^ 4A _ 2), A = 1, 2, ■ ■ • , 

where A is Laplacian on 5"^. They are normalized as follows 



(A.2) 



where 



We consider the following metric (|3.1j ) on 5"^ 

ds^ = cos^ edVt^ + de"^ + sin^ 6'#^ (A.4) 

where 9 G [0, 7r/2] and (f) G [0, 2ti]. 

The SO{A) invariant scalar spherical harmonics of SO{Q) are given by 

y^''=(^,(/>)=CA,fcy^''(^)e^'=^ (A.5) 

where A is a non-negative integer, k G (— A, — A + 2, ■ • • , A — 2, A), ca,A; is a normalization 
constant and y^'^ is given by a hypergeometric function 

y^''^ = sinl^l e 2F1 (^-^(A - \k\), 2 + 1(A + 1 + sin^ . (A.6) 

For the computations in this paper, we need the following scalar spherical harmonics 

i^^'^^ = ^ sin^ ^e±^^^ (A.7) 

y3,±i _ ^ ^2 gin2 e - l)e±*'^ 

4 
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and vector spherical harmonic 

= ^ sin^ 9d(j). (A.8) 
v2 

In order to write down the exphcit form of the dimension A = 2, 3 5*0(4) invari- 
ant chiral primary operators ( p.29|) we write the spherical harmonics in (|A.8|) using the 
following embedding coordinates 



X^ = cos^e, z = l,---,4 

= sin0cos(/), (A.9) 

where Y,t=i ®i = 1- 

Evaluating the harmonics at 6* = 7r/2 we get that y^''^(6' = 7r/2, (p) = CA,ke'^^^, where 
(72,0 = 1/73, C2,±2 = l/2, C3,±i = v/3/4, C3,±3 = l/2v^ (A.IO) 
are obtained from the normalization condition ( [A .21) . 



Appendix B. Supersymmetry of Wilson Loops 

We review here the supersymmetry calculation for the Wilson loop ( |2.38| ), as done in 
2^ and show that they are compatible with those preserved by the surface operator O^,. 



We follow the notations in the Appendix of 0, using 10- dimensional gamma matrices such 
that the variation of the two components of the gauge field gives terms with the gamma 
matrices 7^ and 7"^ and of the two scalars in $ and $ gives 7^ ± i'y^. In addition the 
variation of the extra scalar (p^ gives 7*. 

The supersymmetry variation of the Wilson loop VFeo,i/'o ( P-38| ) gives an expression 
proportional to 

^ sin ■!/;(— 27^ + sin 6*0 sini/^oT^ ~ sin 6*0 cos'i/'o7^) 

+ (cos tjj{i'y^ + sin Oq cos tpo"/^ + sin 6q sin tpol^) + cos 6*07^ j (ei + cos tp'y'^e2 + sin tp'y^e2) ■ 

(B.l) 

Here ei are the parameters of the Poincare supersymmetry transformations and €2 are the 
parameters of the of the conformal supersymmetry transformations. The super symmetries 
preserved by the loop are those where the above expression ( |B.1| ) vanishes for all t/j. One 
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can therefore expand (|B.1|) in Fourier modes and get different equations for each Fourier 
component: 1, sin?/', cos?/', sin2'i/' and cos2'i/'. 

It turns out that requiring aU these projectors to annihilate the spinors is a consistent 
set of equations which have 8 independent solution (for — there are 16). Therefore 
such a Wilson loop is quarter BPS. 

These equations in general relate ei and e2, but we can separate them and find the 
following equation that is satisfied by both spinors 

sin 9q [ sin i^o (7^^ + 7^^) + cos i^o (7^^ - 7'°)] e = . (B.2) 

The supersymmetries unbroken by the surface operator in the plane transverse to the 
loop are the solutions to the equations [0 (and see also the next appendix) 

72389 e = -e , (B.3) 

with either ei or £2- It is easy to check that any solution to (|B.3|) automatically solves 
( [B.2| ) . Therefore the presence of the surface operator Oj: does not break any of the super- 
symmetries left unbroken by the Wilson loop (|2.38| ), they share eight supercharges. 

In the special case of the maximally supersymmetric loop, for Oq = 0, the surface 
operator does break half of the supersymmetries of the Wilson loop, so the combined 
system is quarter supersymmetric, i.e. preserving eight supercharges, like the generic case. 

Appendix C. Supersymmetry of Probe D3-Brane 

In this Appendix, we study the supersymmetries left unbroken by a probe D3-brane 
corresponding to a surface operator Oj:. To avoid the complications of Euclidean su- 
persymmetry, we consider Oj: supported on E = R^'^ in R^'^. We take the metric on 
AdS5 X to be 

^2 3 j-l 4 

ds^ + u^{-dxl + dx\ + dxl + dxl) + c^^'o + H ^^"^^ ^ W ^^^^ 

j=l i=0 i=0 

(C.l) 

The Killing spinor of AdS^ x in this coordinate system is given by |]7^ 

r]2+u;^h{(f),,(f))r]i, (C.2) 



57 



where 7^ are the tangent space gamma matrices and 

771 and 772 are constant ten dimensional complex spinors that have positive and negative 
ten dimensional chirality respectively, 

liiVi = -Vi, 711^2 = 112, (C.4) 

which also satisfy 

771 = (1 - i7oi23)ei, ?72 = (1 + «7oi23)e2, (C.5) 

where ei and 62 are ten dimensional Majorana-Weyl fermions with positive and negative 
chirality. At the boundary, ei is identified with a Poincare supersymmetry parameter and 
€2 with a conformal supersymmetry parameter of the dual gauge theory. 

The embedding of the supersymmetric probe D3 brane (|3.8|) is given in this coordinate 
system by 

sinhwo /> , , , 

uj 2 

where X2 + ixs = r e*^ . In this metric and embedding, the induced metric on the D3- brane 
is AdSs X 

ds^ = cosh^ uq H 2 — {-dxl + dx\) + (i^^ . (C.7) 

V cosh Uq J 

The supersymmetries that are preserved by the D3-brane are those that satisfy 

rKe|_D3 = e\Dz, (C.8) 



where e is the Killing spinor of the AdS^ x background given in (|C.2| ), and is the 
K-symmetry projection matrix 

= -iC-^^Bi dox^ dix" d^xP d^x'' V^^p^ . (C.9) 
For our D3 embedding ( |C.6|) 

-n -n sinh^ Uq sinh , , 

l-C-DBIi^K = i^OlujS ^ n F0123 5 sm(0 — 0o)('^Foitj2 — FoiSfii) 

^ ^ (C.IO) 

5 — cos(0 - 0o)(wFoic^3 + Foi2,^)- 
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The gamma matrices appearing in the last two equations are related to the tangent space 
7 matrices by 

T^=u'j^, r^ = -74, r0 = 79. (c.ii) 

UJ 

Using Cdbi = wcosh uq we may write as 



= A — ( (70149 + sinh^ Wo 70123) - sinhwosin((^ - (f)o){'yoiA2 - 70139) 

cosh ?An V 



cosh uo 

- sinh Uo cos(^ - (po) (70143 + 70129^ 



After some algebra we get from ( |C.8|) 



(1 +72389)ei = 0, 
(1 +72389)e2 = 0. 



(C.12) 



(C.13) 



The preserved super symmetries ( |U.13| ) are the same as those preserved by a maximally 
supersymmetric surface operator Oj: in A/" = 4 SYM (see for e.g. 0). 



Appendix D. Gauge Transformation on V 

As mentioned in the main text, the one- form ( [4.7|) is determined up to an exact 
form by solving the equation dV = ^ *x dz. In order to get a manifestly asymptotically 
AdS^ X geometry, we need to add to V (^3) an exact one-form 

V^V + duj. (D.l) 

The exact form is given up to 0{1/R'^) 



u = a 



2 2sin^i? / 4sin2^i?2 

/ x?i sinSa -f x?9 cosSa — 3x?iX; 2 cosSa — 3x; ix?9 sinSaX /l\ 

Si" ^ — ^ — 



(D.2) 



where we have used the coordinates in (|4.20| ). In computing the correlation functions of 



surface operators Oj^ with local operators using the "bubbling" supergravity solution, we 
have used this expression for V. 
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